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kernels  can  be  exactly  integrated.  We  found  that  the  magnitude  of  the  relative  strain  energy  increases  with 
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ABSTRACT 


Quantum  structures  made  of  components  with  at  least  one  dimension  being  at 
nanoscale,  show  great  potential  for  future  optoelectronic  device  applications.  The  elastic 
fields  in  quantum  structures  affect  their  physical  and  mechanical  properties,  and  also  play 
a  significant  role  in  their  fabrication.  Therefore,  it  is  crucial  that  the  induced  elastic  fields 
in  quantum  structures  be  modeled  accurately  and  efficiently. 

In  Chapter  II,  a  rigorous  analysis  on  the  elastic  and  electric  fields  in  2-dimensional 
quantum  wire  (QWR)  structures  is  presented  using  the  novel  Green’s  functions  and 
related  boundary  element  method  (BEM).  The  elastic  and  electric  fields  in  embedded 
QWR  structures  for  both  the  inclusion  and  inhomogeneity  models  are  investigated.  The 
electric  field  distribution  in  polygonal  QWRs  with  different  sides  is  also  studied  and  it  is 
found  that  the  electric  field  in  triangle  and  square  QWRs  can  be  very  different  to  those  in 
polygonal  QWRs  with  sides  larger  than  4. 

In  Chapter  III,  a  bimaterial  BEM  is  developed  for  the  calculation  of  the  strain  energy 
density  and  the  relative  strain  energy  in  free-standing/embedded  QWR  structures.  The 
required  bimaterial  Green’s  functions  are  derived  in  terms  of  the  Stroh  formalism.  The 
boundary  of  the  QWR  is  discretized  with  constant  elements  for  which  the  involved 
Green’s  function  kernels  can  be  exactly  integrated.  We  found  that  the  magnitude  of  the 
relative  strain  energy  increases  with  increasing  depth  of  the  QWR  with  respect  to  the 
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surface  of  the  substrate.  Strain  energy  density  inside  the  QWR  is  also  plotted  to  show  its 
close  relation  to  the  QWR  shape. 

In  Chapter  IV,  an  analytical  method  for  calculating  the  3-dimensional  quantum  dot 
(QD)  induced  elastic  field  in  the  half-space  substrate  is  presented.  The  QD  is  assumed  to 
be  of  any  polyhedral  shape,  and  its  surface  is  approximated  efficiently  in  terms  of  a 
number  of  flat  triangles  so  that  the  Green’s  function  kernels  can  be  integrated  analytically 
over  the  flat  triangles.  Numerical  examples  are  presented  for  cubic,  pyramidal,  truncated 
pyramidal  and  point  QDs  in  half-space  substrate.  The  strain  energy  distribution  on  the 
surface  of  the  substrate  indicates  clearly  the  strong  influence  of  the  QD  shape  and 
position  on  the  induced  strain  energy.  This  long-range  strain  energy  on  the  surface  is  the 
main  source  for  controlling  and  modulating  the  overgrown  QD  pattern  and  size. 

At  last,  a  detailed  theoretical  calculation  of  the  elastic  and  electric  fields  in  and  around 
nitride-based  QDs  which  are  buried  in  anisotropic  half-space  substrate  is  developed. 
Results  are  presented  for  a  single  dot  as  well  as  coupled  dots.  We  consider  in  detail  the 
case  of  AIN  QDs  in  the  shape  of  hexagonal  truncated-pyramids.  The  calculated  strain  and 
piezoelectric  potential  distributions  induced  by  a  single  QD  are  presented.  Large 
piezoelectric  potential  can  be  observed  in  the  structures.  The  results  are  compared  to 
those  of  simplified  model  in  which  the  QD  is  assumed  to  be  a  point.  Very  similar  trends 
are  observed.  Strain  and  piezoelectric  potential  distributions  induced  by  coupled  QDs  are 
also  shown  along  line  scans  or  on  the  surface. 
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CHAPTER  I 


BACKGROUND 


1.1  Introduction  and  Motivation 

Quantum  confined  nanostructures  have  drawn  great  attention  due  to  their  potential 
application  in  the  fabrication  of  a  wide  variety  of  optoelectronic,  microelectronic  and 
fluorescence  devices,  such  as  photovoltaic  cells,  light  emitting  diodes  and  quantum 
semiconductor  lasers. 

1.1.1  Carrier  Confinement 

Quantum  confined  structures  include  quantum  well  (QW),  quantum  wire  (QWR)  and 
quantum  dot  (QD)  (shown  in  Figure  1.1),  of  which,  the  QWR  and  QD  are  the  most 
popular  topics  and  they  are  also  the  center  subjects  of  this  thesis.  A  key  feature  of  these 
structures  is  carrier  confinement  (electrons  or  holes  are  trapped  in  a  region  with  typical 
dimensions  ranging  from  a  few  nanometers  to  several  hundred  nanometers): 

1.  Quantum  wells,  which  confine  electrons  or  holes  in  one  dimension  and  allow 
free  propagation  in  two  dimensions. 

2.  Quantum  wires,  which  confine  electrons  or  holes  in  two  spatial  dimensions  and 
allow  free  propagation  in  the  third. 

3.  Quantum  dots,  which  confine  electrons  or  holes  in  all  three  spatial  dimensions. 
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Quantum  Quantum  Quantum 

Well  Wire  Dot 


Figure  1.1  Schematic  illustrations  for  QW,  QWR  and  QD 


The  capability  to  confine  carriers  is  based  on  a  difference  in  bandgap  between  two 
materials  (shown  in  Fig.  1.1  the  blue  material  and  its  surrounding  material)  that  are 
brought  in  contact  with  each  other.  Fig.  1.2  shows  the  bandgap  difference  between 
different  materials.  Free  electrons  and  holes  preferentially  stay  in  the  material  with  the 
lowest  bandgap,  where  their  potential  energy  is  minimal.  The  idea  of  reducing 
dimensionality  (Fig.  1.1)  is  leading  to  increasing  confinement.  This  reduction  of 
dimensionality  is  directly  reflected  in  the  density  of  states,  as  shown  in  Fig.  1.3.  The 
density  of  states  has  very  different  shapes  as  a  function  of  energy  depending  on  the 
reduction  of  the  dimensionality:  the  density  of  states  in  a  QW  system  has  a  step  function; 
the  density  of  states  in  a  QWR  system  has  an  inverse  energy  relationship.  Of  all  these 
quantum  structures,  QDs  comprise  the  ultimate  in  carrier  confinement,  trapping  electrons 
and  holes  in  all  three  spatial  dimensions.  It  is  this  fundamental  property  that  gives  rise  to 
the  discrete  energy  spectrum  and  the  density  of  states  in  QD  system  is  continuous  and 
independent  of  energy. 

The  strain  field  in  quantum  structures  can  change  the  bandgap,  cause  trapping  of  holes 
and  electrons,  and  also  cause  electrical  fields  in  some  piezoelectric  materials.  In  short, 


2 


strain  is  a  crucial  factor  to  impact  the  electronic  and  hence  the  optoelectronic  properties 


in  quantum  based  devices. 


Lattice  constant  (nm) 

Figure  1.2  Bandgap  energy  and  lattice  constant  for  different  materials 


quantum  dot 


Figure  1.3  Density  of  states  for  bulk  material,  QW,  QWR  and  QD. 


1.1.2  Fabrication  of  Quantum  Structures 

In  general,  quantum  structures  are  fabricated  with  an  intrinsic  elastic  field  arising  from 
the  lattice  mismatch  between  the  quantum  island  and  substrate  material.  Take  the  InAs 
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quantum  dots  on  GaAs  substrate  as  an  example.  The  lattice  constant  for  InAs  is 

o  o 

«inAs=6.058  A,  while  GaAs  is  acaAs=5.653  A  (Figure  1.4).  Their  natural  lattice  constants 
difference  is  the  so-called  misfit  strain  which  can  be  calculated  as 

y-  _  aGaAs — ajn^L_  g  =y^§  Therefore,  for  InAs/GaAs  quantum  system,  the  misfit  strain 

J  n  J  J 

ulnAs 

is  around  7%. 

When  a  material  is  deposited  on  another  material  (substrate)  with  different  bandgap 
and  lattice  constant,  the  difference  in  lattice  constant  induces  strain  between  the  two 
materials.  Spontaneous  strain  relaxation  lowers  the  potential  energy  of  this  system,  and 
leads  to  the  formation  of  islands  on  the  substrate  (Fig.  1.5).  QDs/QWRs  grown  by  this 
process  are  usually  called  ‘self-organized’  QDs/QWRs.  In  the  growth  process,  the  strain 
plays  a  key  role  in  determining  the  size  uniformity  and  the  spatial  ordering  of  the 
deposited  material. 


/ 
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XX 
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Figure  1.4  Lattice  constants  for  InAs  and  GaAs 
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Self-organization 


Figure  1.5  Schematic  representation  of  quantum  dots  growth. 


In  this  thesis,  we  have  used  two  typical  crystal  lattice  structures  (Fig.  1.6)  which  are 
representative  for  many  quantum  structures:  The  cubic  structure  and  the  hexagonal 
structure.  The  latter  has  a  huge  built-in  electric  field. 


Figure  1.6  Two  typical  crystal  structures:  (a)  cubic  and  (b)  hexagonal 


1.2  Literature  Review 

It  is  well  known  that  the  elastic  and  electric  fields  in  quantum  structures  substantially 
modify  the  electronic  band  structure  and  thus  strongly  affect  the  performance  of  the 
electronic  devices.  One  could  control/optimize  the  electronic  and  optical  properties  of 
quantum  structures  by  altering  the  elastic  and  electric  fields  [ConOO,  Med02,  Not95]. 
Also,  a  detailed  knowledge  of  the  strain  field  can  serve  as  a  useful  tool  to  understand  the 
growth  of  these  structures  in  order  to  control  the  size,  shape  and  distribution  of  the 
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wires/dots,  to  achieve  good  repeatability.  Hence,  it  is  crucial  that  the  induced  elastic  and 
electric  fields  should  be  modeled  accurately  and  efficiently  in  order  to  obtain  a  well 
ordered  quantum  structure  and  to  improve  the  performance  of  the  electronic  devices. 

Mathematical  and  numerical  simulations  play  an  important  role  in  the  study  of 
quantum  confined  nanostructures.  At  the  research  point  of  view  for  exploring  the  physics, 
computational  simulation  is  an  instrument  in  interpreting  related  data  and  guiding  for 
future  developments.  At  the  engineering  stage  of  optimization,  it  can  provide  quantitative 
analysis  quickly  and  economically. 

There  are  a  variety  of  techniques  been  developed  by  many  researchers  to  study  the 
elastic  and  electric  fields  in  the  quantum  nanostructures,  such  as  the  finite  element 
method  (FEM)  and  finite  difference  method  [Ben96,  Kre99,  Pei03,  Gru95,  Liu02,  Jon06], 
atomistic  approach  [Pry98,  CalOO,  KikOl,  MakOl],  continuum  approaches  [JogOO,  GlaOl, 
Zha05].  Jiang  and  Singh  determined  the  strain  field  for  InAs  QD  embedded  in  GaAs 
substrate  using  a  valence  force  field  model  (one  of  the  atomistic  models)  [Jia98].  But 
when  a  large  number  of  atoms  are  required  for  accuracy  in  the  study,  this  method  appears 
computationally  intensive  and  subjects  to  convergence  problems.  Furthermore,  the 
atomistic  approach  depends  crucially  on  the  quality  of  the  interatomic  potential  model. 
The  key  advantage  for  this  technique  is  that  displacements  are  evaluated  at  the  atomistic 
level  which  can  provide  potentially  greater  accuracy  than  other  methods  at  the  QD  edges 
and  the  QD/substrate  interface. 

The  FEM  approach  is  the  most  commonly  used  method  in  engineering  fields. 
Grundmann  et  al  used  the  FEM  to  analyze  a  pyramidal  dot  [Gru95].  Pei  et  al.  investigated 
the  elastic  fields  in  the  self-organized  quantum  dot  structures  by  a  three-dimensional 
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FEM  [Pei03].  Liu  et  al.,  also  using  the  FEM,  reported  the  stress  and  strain  of  the 
heterosystem  of  InAs  quantum  dots  embedded  in  GaAs  [Liu02].  Jonsdottir  et  al  used  the 
FEM  to  assess  the  effect  of  lattice  mismatch,  dot  volume  and  dot  surface  contact  area  on 
the  induced  stresses  and  elastic  energies  [Jon06].  The  FEM  approach  could  be  also 
computationally  intensive  since  it  requires  discretization  of  the  problem  domain. 

Analytical  continuum  methods  for  calculating  the  strain  distribution  induced  by  a 
uniformly  lattice  mismatched  inclusion  were  based  on  the  classic  work  of  Eshelby 
[Esh57].  Elastic  field  induced  by  QDs  are  determined  by  performing  an  integral  over  the 
surface  of  QDs.  Davies  extended  the  elastic  field  evaluation  to  electroelastic  field  and 
determined  the  piezoelectric  potential  [Dav98].  The  advantages  of  the  analytical 
continuum  methods  are  quick  and  portable,  and  accurate.  However,  there  are  some 
restrictions  associated  with  them.  These  include  assumption  of  isotropic  material 
properties,  same  properties  for  both  QD  and  substrate,  and  of  the  infinite  substrate. 

It  is  well-known  that  any  pure  continuum  method  cannot  account  for  atomistic  scale 
distortion.  Andreev  et  al.  then  proposed  an  exact  analytical  solution  for  the  lattice 
relaxation  associated  with  a  variation  of  lattice  constant  [And99].  The  analytical  solutions, 
in  particular  those  related  to  the  Green’s  function  solutions,  have  been  developed  and 
applied  to  the  quantum  structure  modeling  [PanOl].  Because  of  their  accuracy  and 
efficiency,  these  analytical  methods,  particularly  the  Green’s  function  method,  could  be 
more  appealing  to  experimentalists  and  physical  device  designers  in  the  studies  of 
quantum  structures  [PeaOO].  More  recently,  the  multiscale  Green’s  functions  have  been 
also  introduced  to  model  QD  induced  strain  field  [Yan05,  Rea07],  which  has  been  shown 
to  be  very  computationally  efficient. 
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From  the  above  literature  review,  it  can  be  concluded  that  at  least  one  or  more  of  the 


following  simplifications  are  made  through  the  calculations. 

1.  Purely  elastic  field.  Most  previous  studies  considered  only  the  strain  field  in  the 
quantum  structures  although  the  electric  field  is  also  a  crucial  factor. 

2.  An  infinite  substrate.  Most  previous  analyses  assumed  that  the  QD  was  buried  in 
an  infinite  matrix  although  the  QD  could  be  buried  immediately  below  a  free 
surface. 

3.  A  single  isolated  quantum  wire  or  dot. 

4.  A  simple  shaped  quantum  wire  or  dot. 

Recently,  Pan  [Pan02c]  has  employed  the  equivalent  body-force  concept  of  Eshelby 
inclusion  problem  [Esh57,  Mur87]  and  succeeded  in  extending  the  Green’s  function 
approach  to  anisotropic  substrates  [Jog03,  PanOl,  Pan02a,  Pan02b,  Yan03].  The  new 
Green’s  function  takes  into  account  quantum  wire  or  dot  of  arbitrary  shape  and  substrate 
with  the  most  general  class  of  anisotropy  and  piezoelectricity.  In  this  thesis,  we  use  the 
extended  Green’s  function  and  the  corresponding  BEM  to  study  the  elastic  and  electric 
fields  in  or  around  quantum  wire  or  dot  systems. 

1.3  Thesis  Outline 

The  thesis  is  organized  as  follows:  In  Chapter  II  we  propose  a  simple  BEM 
formulation  to  investigate  the  elastic  and  electric  fields  in  QWR  structures.  In  Chapter  III 
the  bimaterial  Green’s  functions  in  anisotropic  elastic  media  are  derived  in  terms  of  the 
Stroll  formalism.  The  strain  energy  density  and  the  relative  strain  energy  in  free-standing 
QWR  or  in  QWR  embedded  in  half-space  substrate  are  calculated.  An  analytical  method 
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for  calculating  the  QD-induced  strain  fields  in  half-space  substrates  are  presented  in 
Chapter  IV.  In  Chapter  V,  a  detailed  calculation  of  the  elastic  and  electric  fields  in  and 
around  nitride -based  QDs  is  developed.  Finally,  the  conclusions  and  contributions  of  this 
research,  as  well  as  future  research  opportunities  are  discussed  in  Chapter  VI. 
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CHAPTER  II 


ELASTIC  AND  ELECTRIC  FIELDS  INDUCED  BY  QWR  IN  FULL-PLANE: 
INCLUSION  VS.  INHOMOGENEITY 


2.1  Summary 

The  well-known  Eshelby  inclusion  method  [Esh57]  has  been  successfully  applied  to 
study  the  induced  elastic  and  electric  fields  in  QWR  structures.  In  this  treatment  the 
lattice  misfit  between  the  QWR  and  substrate  is  allowed  but  the  QWR  material  is 
assumed  to  be  the  same  as  its  substrate  [Fau96,  Dav98,  RuOO,  Gla03,  Pan04,  Mar06].  The 
most  advantage  of  the  Eshelby  inclusion  method  is  its  simplicity  by  which  the  induced 
fields  can  be  estimated  analytically. 

Recently,  the  accuracy  of  the  inclusion  model  has  been  studied  using  the 
inhomogeneity  model.  Different  material  property  for  the  QWR  and  substrate  was 
applied  [Yan02,  Pan05],  instead  of  assuming  the  same  material  property.  For  instance, 
the  issue  of  the  homogeneous  inclusion  vs.  inhomogeneity  in  the  context  of  QWR 
structures  was  studied  in  detail  for  the  purely  elastic  structures  by  Pan  et  al  [Pan05]. 
Numerical  examples  were  given  for  InAs/GaAs  QWR  structures  in  both  (001)  and  (111) 
growth  directions  for  QWRs  in  shape  of  square  and  trapezoid.  There  are  some  main 
observations  from  their  analysis:  (i)  In  the  substrate  and  far  away  from  the  QWR,  both 
the  inclusion  and  inhomogeneity  models  give  similar  results,  (ii)  For  points  within  or  near 
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the  QWR,  the  difference  between  the  inclusion  and  inhomogeneity  models  can  be  as  high 
as  10%  for  these  materials,  (iii)  Although  the  singular  behavior  near  the  comers  of  the 
QWR  looks  similar  for  both  the  inclusion  and  inhomogeneity  models,  the  amplitudes  of 
the  singularity  are  different.  So  far,  however,  the  difference  on  the  electric  field  based  on 
both  the  inclusion  and  inhomogeneity  in  the  corresponding  InAs/GaAs  (111)  QWRs  and 
in  QWR  made  of  any  nitride  III  group  material  has  not  been  reported  in  any  literature, 
which  is  the  main  motivation  of  this  study. 

This  chapter  is  organized  as  follows.  In  Section  2.2,  we  describe  the  associated  basic 
equations.  In  Section  2.3,  we  present  the  problem  to  be  solved.  In  Section  2.4,  the  BEM 
and  the  corresponding  constant  element  discretization  are  presented.  While  in  Section  2.5 
we  present  various  numerical  examples,  conclusions  are  drawn  in  Section  2.6. 

2.2  Basic  Equations 

The  basic  equations  for  a  linear  and  generally  anisotropic  piezoelectric  solid  under  a 
static  deformation  condition  conclude  the  constitutive  relations 

&ij  ~  CijmnYmn  ~ekij^k  (2- la) 

~  eijk  Y jk  ^ijEj  (2.  lb) 

where  al}  is  stress  and  D\  is  electric  displacement; is  strain  and  E\  is  electric  field;  Cijmn, 
and  are  elastic  moduli,  the  piezoelectric  constants  and  dielectric  constants 
respectively. 

The  strain-displacement  relations 
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(2.2) 


n=\(uu+uu) 

E,  =-</>,  (2.3) 

where  w,  and  0  are  the  elastic  displacement  and  electric  potential  respectively.  And 
the  equilibrium  equations 

<r„+fj=0  (2-4) 


-q  =  0 


where /;  and  q  are  body  force  and  electric  charge  respectively. 


For  the  seek  of  convenience,  we  adopt  the  extended  notations  for  the  stress  and 
electric  displacement,  elastic  strain  and  electric  field,  elastic  displacement  and  electric 
potential,  body  force  and  electric  charge,  and  all  the  three  kinds  of  material  constants  in 
piezoelectric  problems  [Bar75,  Dun93]  as 


n  = 

Uj  = 

F,= 


K’ 

J  =  j  =  1,2,3 

l  A, 

J  =  4 

\r„ . 

/  =  1,2,  3 

,  1  =  4 

I  =  i  =  1,2,3 

. 

1=4 

[fn 

I  =  i  =  1,2,3 

\-q , 

1  =  4 

(2.6) 


(2.7) 


(2.8) 


(2.9) 
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(2.10) 


Cm ,  J,K  =  j,k  =  1,2,3 
e„y,  7  =  7  =  1, 2, 3;  K  =  4 

I  =  < 

V-  :  r  1^1 

em,  7  =  4;  K  =  k  =  1,2,3 
-en,  J  =  K  =  4 

In  this  definition,  it  can  be  noted  that  the  original  symbols  are  kept  just  with  the  fourth 
components  related  to  electric  fields.  In  terms  of  this  extended  notation,  the  constitutive 
relations  can  be  simplified  to  a  single  equation 

O a  =CiJKlrKi  (2.11) 

Similarly,  the  equilibrium  equations  can  be  unified  as 

<rUj+Fj=  0  (2-12) 


2.3  Problem  Description  for  QWR  in  Full-plane 

Consider  a  general  QWR  problem,  which  is  illustrated  in  Figure  2.1.  The  QWR  can  be 
arbitrarily  shaped.  Let  us  assume  that  the  general  misfit  strain  y*tj  (1=1 ,2,3,4  and  j=l ,2,3) 

is  uniform  within  the  QWR  domain  and  is  zero  outside. 

The  total  extended  strain  is  the  sum  of 

n  =  n+rl  (2.13) 


where  y^  the  extended  strain  that  appears  in  the  generalized  constitutive  relation  in  Eq. 
(2.11)  as 


CukiYkI 


(2.14a) 


or 


—  CiJKl  ( Yki  XYki) 


(2.14b) 


13 


In  Eq.  (2.14),  %  is  equal  to  1  if  the  observation  point  is  within  the  QWR  domain  V  and 


to  0  outside. 


z 


Figure  2.1  An  arbitrarily  shaped  polygon  QWR  inclusion/inhomogeneity  with 
eigenstrain  yjj  in  an  anisotropic  substrate. 

The  interface  between  the  QWR  and  matrix  is  labeled  s.  We  also  denote  C"JKI  and 

C'"JKl  as  the  general  moduli  of  the  QWR  and  matrix  materials,  respectively  (Figure  2.1). 

For  the  homogeneous  inclusion  problem,  CwUKl  =C’’jKl . 

Assuming  that  there  is  no  body  force  and  electric  charge  within  the  QWR  system,  one 
can  easily  show  that,  for  the  matrix  domain, 


Cm  m  _  r\ 

uvMv  U 


(2.15) 


iJKluK,li 


and  for  the  QWR  domain, 


(2.16) 


It  is  clear  that  the  right  hand  side  of  Eq.  (2.16)  is  equivalent  to  a  body  force  defined  as 


(2.17) 
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which  is  also  called  the  equivalent  body  force  of  the  eigenstrain  [Mur87,  Pan04a].  This 
equivalent  body  force  will  be  employed  in  the  next  section  to  convert  the  contribution  of 
the  eigenstrain  to  a  boundary  integral  along  the  interface  of  the  QWR  and  its  substrate. 
Again,  the  superscripts  m  and  w  denote  quantities  associated  with  the  matrix  and  QWR, 
respectively. 

2.4  Boundary  Integral  Equations  in  Terms  of  Green’s  Function 

To  solve  the  problem  in  Figure  2.1,  BEM  is  applied  to  both  the  QWR  and  its 
matrix/substrate.  The  boundary  integral  formulation  can  be  expressed  as  [Yan02,  Pan05] 

K  (X)u';  (x) = \s  [t /;;  (x,x>;  (x)  -  t;;  (x,  x)«;  (x)^(x)  (2. 1 8) 

for  the  matrix,  and 

bu  (x)«;(x) = \s  [uz  (x,x)fr;(x)+ /;(x))-  r;(x,x)M;(x)^(x)  (2. 19) 

for  the  QWR. 

In  Eqs.  (2.18)  and  (2.19),  jc(.v,  z)  and  X(X,  Z)  are  the  coordinates  of  the  field  and  source 

points.  Coefficients  bu  is  equal  to  8u  if  X  is  an  interior  point  and  du  at  a  smooth 

boundary  point.  For  points  at  complicated  geometry  locations,  these  coefficients  can  be 
determined  by  the  rigid-body  motion  method  [Bre92].  tj  and  u\  are  the  traction  and 
displacement  components,  respectively.  Furthermore,  in  Eq.  (2.19),  fj  is  the  traction 
induced  by  the  misfit  eigenstrain  inside  the  QWR  whose  expression  is  given  by  Eq. 
(2.17). 
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The  Green’s  functions  Uu  and  T\\  in  Eqs.  (2.18)  and  (2.19)  are  taken  to  be  the  special 


2D  Green’s  functions  for  the  full  plane  [Tin96,  Pan02a].  The  indexes  7  and  J  indicate  the 
/-th  Green’s  (general)  displacement/traction  (at  x)  in  response  to  a  (general)  line-force  in 
the  7-th  direction  (applied  at  X).  Note  that  the  Green’s  functions  are  in  exact  closed  form, 
and  thus  their  integration  over  constant  elements  can  be  carried  out  exactly  as  discussed 
below.  This  is  computationally  desirable  as  it  is  very  efficient  and  accurate  for  the 
calculation. 


z 


Figure  2.2  QWR  boundary  discretization 


To  employ  the  above  Green’s  Functions,  we  divide  the  boundary  of  QWR  into  N 
segments  with  the  n-th  element  being  labeled  as  sn  (shown  in  figure  2.2).  The  constant 
values  Mjn  and  ?jn  on  the  «-th  element  equal  to  the  displacement  and  traction  at  the  center 
of  the  element.  Under  this  assumption,  the  boundary  integral  equations  (2.18)  and  (2.19) 
for  the  surrounding  matrix  and  QWR  domains  are  reduced  to  the  following  algebraic 
equations 


N 


N 


(2.20) 
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for  the  matrix,  and 


*„«;  +E<[  =Z<|,  <2.21) 

n= 1  ”  «=1  ” 

for  the  QWR. 

It  is  obvious  that  the  difference  between  Eqs.  (2.20)  and  (2.21)  is  the  traction  induced 
by  the  misfit  eigenstrain  inside  the  QWR  in  Eq.  (2.21). 

The  problem  now  is  reduced  to  find  the  suitable  Green’s  functions  U\j  and  Tu,  as  well 
as  their  integrals  over  each  element  ,vn,  which  are  the  kernel  functions  in  these  equations. 
We  now  present  the  analytical  integration  of  these  Green’s  functions  over  an  arbitrary 
constant  element. 

To  carry  out  the  line  integration  of  the  Green’s  functions  over  a  constant  element,  we 
first  look  at  the  Green’s  functions  appeared  in  Eqs.  (2.20)  and  (2.21).  They  can  be  simply 
expressed  as  [Pan04b,  Tin96,  Pan02b] 

Uu  (x,X)  =  -Im{Ays  ln(zs  -  sR)A/R}  (2.22) 

n 

Tu  (x,  X)  =  -  Im j Bjr  PRfll~n 3  Am  }  (2.23) 

7T  y  zR  Sr  J 

where  “Im”  stands  for  the  imaginary  part  of  the  complex  value,  Au  and  Bu  are  two 
constant  matrices  related  to  the  material  property  only  [Tin96],  n\  and  773  are  the  unit 
outward  normal  along  x-  and  z-directions,  pR  (R=l,2,3,4)  are  the  Stroh  eigenvalues,  and 
ZR=x+pRz  and  s R=X+pRZ  are  related  to  the  field  x(.v,  z)  and  source  X(X.  Z)  points, 
respectively. 
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The  Stroll  eigenvalues  and  eigenmatrices  involved  in  Eqs.  (2.22)  and  (2.23)  are 
obtained  by  solving  the  following  equations.  The  eigenvalue  p  and  the  corresponding 
eigenvector  a  are  solved  from  the  equation 

[Q  +  p(R  +  RT)  +  p2T]a=  0  (2.24) 

where  the  superscript  T  denotes  matrix  transpose.  And 

Gu=Cmi,  fiu=Ciu3,  Tu=C3m  (2.25) 

where  CL|Ki  are  the  elastic  and  electric  moduli  defined  in  Eq.  (2.10). 

After  we  get  eigenvector  a,  the  eigenvector  b  can  be  obtained  from 

b  =  (RT  +PT)a  =  --(Q  +  PR)a  (2.26) 

P 

Denoting  by  pk,  a k  and  bk  (k=l,2,3,...,8)  the  eigenvalues  and  associated  eigenvectors, 
we  order  them  in  such  a  way  so  that 

lmpi>0,  pI+4  =  pr ,  aI+4  =a, ,  bI+4  =b,  (1=1, 2, 3, 4)  (2.27) 

where  an  overbar  denotes  the  complex  conjugate,  and  the  matrices  A  and  B  are 

A=[ai,a2,a3,a4],  B=[bi,b2,fa,b4\  (2.28) 

We  assume  that  p\  are  distinct  eigenvalues,  a\  and  b\  satisfy  the  normalization  relation 

bTIaJ+aTIbJ=8IJ  (2.29) 

with  Su  being  the  4  by  4  identity  matrix.  Remarked  that  repeated  eigenvalues  p\  can  be 
avoided  by  using  slightly  perturbed  material  coefficients  with  negligible  errors. 
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We  now  consider  the  integrals  of  Green’s  functions  over  a  constant  element.  Define  a 
generic  line  segment  for  any  constant  element  along  the  interface  in  the  {x,z)~ plane, 

starting  from  point  1  (jci,  z\)  and  ending  at  point  2  (x2,  Zi).  In  terms  of  the  parameter  t 

(0<t<l),  any  constant  line  element  can  be  expressed  as  [Pan04b] 

x  =  X\  +  (x2  -xx)t  (2.30a) 

z  =  z,  +(z2  -Z\)t  (2.30b) 

The  outward  normal  components  along  the  line  segment  are  constant,  given  by 

n1=(z2-z1)/l;  n3  =  ~(x2  - xx)  1 1  (2.31) 

where  /=  2J(x2  -  xx)2  +  (z2  ~  zx)2  is  the  length  of  the  line  segment.  The  variable  of 
integration  ds  in  Eqs.  (2.20)  and  (2.21)  is  as  ds=/df. 

Consider  the  corresponding  integrals  J  U uds  and  j  Tuds  in  Eqs.  (2.20)  and  (2.21). 

Substitute  the  Eqs.  (2.22)  and  (2.23)  into  the  above  two  expressions  to  get 

J  Uuds  =  /J  (—  Im{A/s  ln(zA>  - sR)AIR})dt  (2.32) 

f  Tuds  =  /  f  (—  1m  |  B,r  P«n'~ni  Aik  \)dt  (2.33) 

Js"  Js"  ^  l  zR-sR  J 

Since  pR,  AJR  and  B}R  are  constant  in  each  segment,  Eqs.  (2.32)  and  (2.33)  consist  of 
two  integrals  involving  the  following  integration.  For  the  first  integral,  we  define  it,  being 
a  function  of  the  source  point  X  =(X,  Z),  as 

hR(X,Z)=f]n(zRsR)dt  (2.34) 

JO 

or, 
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hR(X,Z)  =  J  ln{[(x2  -Xy)  +  pR(z2-zl)]t  +  [(x1  +  Pgzj-Sgftdt  (2.35) 


Integration  of  this  expression  gives 


hR(X,Z)  = 


_  (*1  +  PRzl)-sh 


(x2-x1)  +  pR(z2-zl) 


-In 


x2  +  PrZ2~sr 
X\  +  PrZ,-Sr 


+  ln [x2  + pRz2~ sR]-l  (2.36) 


Similarly,  we  define  the  second  integral  as 


gR(X,Z)sj‘ — — 

JO  y  —  o 


(2.37) 


Integration  of  this  expression  gives 


gR(X,Z)- 


1 


(x2-x1)  +  pR(z2-zl) 


-In 


x2  +  pRZ2~sR 

X1  +  PrZi-Sr 


(2.38) 


Finally,  the  integral  of  the  displacement  and  traction  Green’s  function  from  the 
contribution  of  a  constant  boundary  element  can  be  obtained  in  the  closed  form  as 

f  UIJds  =  -lm[AJRhR(X,Z)AIR\  (2.39) 

Jsn  n 

f  TlJds  =  -lm[BJRgR(X,Z)AIR]  (2.40) 

Js»  n 

Therefore,  using  the  constant  element  discretization  the  two  boundary  integral 
equations  (2.18)  and  (2.19)  for  the  QWR  and  matrix/substrate  can  be  cast  into  a  system  of 
algebraic  equations  for  the  interface  points.  In  matrix  form,  they  can  be  expressed  as 


mw  _(T"  +bi)uw  =r 

(2.41) 

U"'t"'-(T"' +M)u"'  =0 

(2.42) 
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where  the  coefficient  matrices  U  and  T  are  the  exact  integrals  of  Green’s  functions  on 
each  constant  element  given  in  Eqs.  (2.39)  and  (2.40),  I  is  the  identity  matrix  with 
I(i,j)=«5j,  and  u  and  t  are  the  general  displacement  and  traction  vectors  in  the  center  of 
each  constant  element.  The  right-hand  side  term/"'  in  Eq.  (2.41)  is  the  general  equivalent 
force  corresponding  to  the  misfit  eigenstrain  within  the  QWR. 


We  assume  that  the  matrix  and  QWR  are  perfectly  bonded  along  the  interface  s ,  that  is, 
the  continuity  conditions  um  =  ww  and  /'"=  -  f  hold  on  the  interface.  We  already  know  that 

bu=—  du  when  source  point  X  is  on  boundary.  Then  Eqs.  (2.20)  and  (2.21)  can  be 


expressed  in  the  matrix  form 


Tm  +  — 
2 

T"  +- 
2 


4 

4 


Um 
—U " 


r  Mwi 

“  0  “ 

u  \ 

r_ 

(2.43) 


In  Eq.  (2.43),  the  number  of  unknowns  is  identical  to  the  number  of  equations  so  all 
the  nodal  (general)  displacements  and  tractions  can  be  determined  by  Gauss  elimination 
method.  Furthermore,  making  use  of  the  general  Somigliana’s  identity,  the  general 
displacement  at  any  location  within  the  QWR  can  be  easily  obtained  using  the  solved 
interface  values 


=  E<1,  t^nc  +  fyf  y,"dr)/;-£([ 

n= 1  n  n= 1  ”  n= 1  n 

while  the  displacement  at  any  location  within  matrix  can  be  obtained  as 

V  N 

b„u,  =£(J,  u-dTK-^l  T-'drK 

n= 1  n  n= 1  n 


(2.44) 


(2.45) 
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In  order  to  find  the  strain  and  electric  fields,  we  need  to  take  the  derivation  of 
displacement  with  respect  to  the  field  point  x  (refer  to  Eqs.  (2.2)  and  (2.3)).  We  first  write 
Eqs.  (2.22)  and  (2.23)  as 

UIJ(x,X)  =  -lm{AJRdR(x,X)AIR}  (2.46) 

n 

TIJ(x,X)  =  —Im{BJR(pRnl -n,)eR(x,X)AIR}  (2.47) 

71 

where 

dR(x,X)  =  \n(zR-sR),  eR(x,X)  =  — - —  (2.48) 

Zr~sr 

The  strain  field  can  be  obtained  by  taking  the  derivative  of  Eqs.  (2.46)  and  (2.47).  Take  a 
point  in  matrix  as  an  example 


1  V  N  N  N 

T„  =  -  { [£(£  u-jry;,  -  £(|  T‘„dT)ul  ] + [£(£  -  £  ([ 

(2.49) 


n= 1 


n= 1 


«=1 


71=1 


W  Al 

71=1  ”  71=1  " 


(2.50) 


where  p,g=l  and  3,  and 


( *.  * )  =  7  to  [AmdRM  (x,  X ) A,„ }  (2,51) 

TC 

T,j,P(qM,X)  =  -lm{BJR{pRnx - n3)eRp(q)(x,X)AIR }  (2.52) 

7T 


where 


dRl(x,X)  = 


z»  ~  sK 


dR3(x,X)  = 


Pr 


ZR  -  s„ 


(2.53) 
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1 


(2.54) 


eRl(x,X)  =  -- -  e  (Xtx)  =  -~  Pr 

(zr~sr)  (zr-sr) 


Furthermore,  utilizing  the  basic  Eq.  (2.1),  all  the  internal  elastic  and  electric  response 
in  the  matrix  and  QWR  can  be  calculated. 

In  summary,  we  have  derived  the  exact  boundary  integral  equations  for  the  QWR  and 
matrix  domains  by  assuming  constant  elements  along  their  interface.  These  equations  can 
be  used  to  find  the  elastic  and  piezoelectric  responses  along  the  interface  and  at  any 
location  within  the  QWR  and  its  surrounding  matrix.  Applications  of  these  solutions  to 
QWR  systems  are  discussed  in  the  next  section. 


2.5  Numerical  Examples 

Materials  used  in  this  numerical  section  consist  of  two  QWR  systems:  One  is 
InAs/GaAs  with  (001)-  and  (11  Indirections,  and  the  other  InN/AIN  with  (1000)-  and 
(000 indirections.  Both  the  inclusion  and  inhomogeneity  models  are  studied.  For  the 
corresponding  inclusion  model,  the  QWR  materials  will  be  GaAs,  instead  of  InAs,  in  the 
first  system,  and  AIN,  instead  of  InN,  in  the  second  system.  The  material  properties  and 
eigenstrains  within  the  QWR  in  the  material  coordinates  are  listed  in  Tables  1  and  2 
[Pan02a,  Pan02b] . 

For  InAs  (111)  and  GaAs  (111),  the  coordinate  x-axis  is  along  [11-2],  y-axis  along  [- 
110],  and  z-axis  along  [111]  directions  of  the  crystalline  [Pan02b].  Their  material 
properties  can  be  obtained  by  simple  coordinate  rotations.  Similarly,  for  InN  (1000)  and 
AIN  (1000),  the  material  properties  can  be  obtained  by  simply  exchanging  the  coordinate 
directions. 
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Table  2.1  Material  properties  &  misfit  strains  in  InAs(001)/GaAs  (001) 


InAs 

GaAs 

CU=C22=  C33  (GPa) 

83.29 

118.8 

c12=c13=  c23 

45.26 

53.8 

C44=C55=  C66 

39.59 

59.4 

Cl4=  625=  636  (C/m“) 

-0.0456 

-0.16 

£u=£22=£33  (10  9C2/Nm2) 

0.1345808 

0.110675 

Yu  =  Y 22~  Y33 

0.07 

Table  2.2  Material  properties  &  misfit  strains  in  InN(0001)/AlN  (0001) 


InN 

AIN 

C11=C22  (GPa) 

223.0 

396.0 

C33 

224.0 

373.0 

Ci2 

115.0 

137.0 

Ci3=C23 

92.0 

108.0 

C44=C55 

48.0 

116.0 

C66 

54.0 

129.5 

ei5=e24  (C/m2) 

-0.22 

-0.48 

e3i=e32 

-0.57 

-0.58 

633 

0.97 

1.55 

£,,=£22  (10"12C2/Nm2) 

132.81 

79.686 

£33 

132.81 

97.372 

*  * 

Yu  =  Y22 

0.1357 

* 

Yu 

0.1267 

2.5.1  Circular  and  Elliptic  QWRs 

The  first  example  is  for  circular  and  elliptic  QWRs  in  an  infinite  substrate  (Figure  2.3). 
For  this  case,  exact  solutions  can  be  obtained  for  both  the  inclusion  and  inhomogeneity 
models. 
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Figure  2.3  A  circle  QWR  of  radius  r=10nm  in  an  infinite  substrate  (a),  and  an  ellipse 
QWR  with  length  of  the  semi-major  axis  a=20nm  in  horizontal  v-direction  and  the  semi¬ 
minor  axis  /?=5nm  in  vertical  z-dircction  in  an  infinite  substrate  in  (b). 

Tables  2. 3-2. 6  compare  the  strain  and  electric  fields  inside  both  the  circular  and 
elliptic  inclusion  and  inhomogeneity  (Figure  2.3),  using  the  material  properties  and 
eigenstrain  listed  in  Tables  2.1  and  2.2.  We  point  out  that  for  these  QWR  shapes,  the 
induced  fields  inside  the  QWR  are  uniform  for  both  inclusion  and  inhomogeneity  models. 
These  results  can  be  obtained  using  the  analytical  solution  for  the  QWR  inclusion 
problem  [Pan04a],  combined  with  the  Eshelby  inhomogeneity  method  [Esh57,  Esh61, 
Mur87].  Furthermore,  we  have  also  used  our  BEM  formulation  presented  above  for  these 
models.  In  doing  so,  we  have  mutually  checked  our  analytical  and  numerical  solutions. 

It  is  observed  from  Tables  2. 3-2. 6  that  the  relative  error,  defined  as  (inclusion  solution 
-  inhomogeneity  solution)/  (inhomogeneity  solution),  for  the  strain  and  electric  field 
(10  V/m)  between  the  inclusion  and  inhomogeneity  models  is  generally  around  10%, 
which  is  consistent  with  the  recent  prediction  for  the  purely  elastic  QWR  semiconductor 
[Yan02,  Pan05].  We  also  note  that  for  the  elliptical  QWR  case,  the  horizontal  strain  yxx 
based  on  the  inhomogeneity  model  can  be  much  different  than  those  based  on  the 
inclusion  model.  Perhaps  the  most  important  feature  is  on  the  difference  of  the  electric 
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fields  in  the  InAs/GaAs  (111)  system.  It  is  observed  that  the  electric  field  difference 
based  on  the  inclusion  and  inhomogeneity  models  can  be  more  than  one  order  of  the 
magnitude  (Table  2.4).  This  special  feature  has  not  been  reported  in  any  previous 
investigation  and  will  be  discussed  again  in  the  following  two  examples. 

Table  2.3  Strains  in  inclusion  GaAs  (001)  and  inhomogeneity  InAs(001)/GaAs  (001) 


Inclusion 

Inhomogeneity 

Relative  Error  (%) 

Circle 

Yxx— Yzz 

0.0612 

0.0566 

9 

Ellipse 

Yxx 

0.0208 

0.0134 

55 

Yzz 

0.1058 

0.1133 

7 

Table  2.4  Strains  in  inclusion  GaAs  (111)  and  inhomogeneity  InAs(l  1  l)/GaAs  (1 1 1) 


Inclusion 

Inhomogeneity 

Relative  Error  (%) 

Circle 

Yxx 

0.0596 

0.0551 

8 

Yzz 

0.0543 

0.0500 

9 

2Yxz 

-0.0149 

-0.0145 

2 

Ex 

(xl07V/m) 

4.020 

0.2738 

1368 

Ez 

(xl07V/m) 

-2.843 

-0.1936 

1368 

Ellipse 

Yxx 

0.0260 

0.0196 

33 

Yzz 

0.0829 

0.0861 

4 

2Yxz 

-0.0081 

-0.0075 

9 

Ex 

(xl07V/m) 

1.235 

0.4773 

176 

Ez 

(xl07V/m) 

-9.770 

2.074 

571 

Table  2.5  Strain  in  inclusion  A1N(000 


)  and  inhomogeneity  InN(0001)/AlN(0001) 


Inclusion 

Inhomogeneity 

Relative  Error  (%) 

Circle 

Yxx 

0.1184 

0.1063 

11 

Yzz 

0.0984 

0.0876 

12 

Ez 

(xl07V/m) 

-47.05 

-42.65 

10 

Ellipse 

Yxx 

0.0515 

0.0310 

66 

Yzz 

0.1555 

0.1746 

11 

Ez 

-159.0 

-130.1 

22 
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Table  2.6  Strain  in  inclusion  AIN(IOOO)  and  inhomogeneity  InN(1000)/AlN(1000) 


Inclusion 

Inhomogeneity 

Relative  Error  (%) 

Yxx 

0.0984 

0.0876 

12 

Circle 

Yzz 

0.1184 

0.1063 

11 

Ex 

(xl07V/m) 

-47.05 

-42.65 

10 

Yxx 

0.0403 

0.0223 

81 

Ellipse 

Yzz 

0.1797 

0.2008 

11 

Ex 

(xl07V/m) 

2.435 

2.780 

12 

2.5.2  Square  QWR  in  GaAs  and  hexagon  QWR  in  AIN 

The  second  example  is  for  a  square  QWR  in  GaAs  and  a  hexagon  QWR  in  AIN 
(Figure  2.4)  where  the  induced  strain  and  electric  fields  are  presented  along  the  horizontal 
and  diagonal  (inclined  in  hexagon)  lines  (dashed  lines  in  Figure  2.4). 


Figure  2.4  A  square  QWR  in  GaAs  (a)  and  a  hexagon  QWR  in  AIN  (b). 


We  first  test  the  convergence  of  our  program.  Take  the  square  InAs/GaAs  (111)  as  an 
example.  The  interface  of  the  square  is  discretized  to  100  or  120  segments.  The  strain 
(Table  2.7)  and  electric  field  (Table  2.8)  at  some  locations  along  horizontal  line  in  Fig. 
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2.4(a)  are  listed.  The  results  show  that  the  errors  are  negligible  when  using  different  mesh 


sizes  and  our  program  has  good  convergence  when  the  mesh  size  reaches  100. 


able  2.7  Strain  yxx  in  inhomogeneity  InAs/GaAs(l  1 1)  for  different  mesh  sizes 


^^^Alesh  size 
x  (nm) 

100 

120 

Relative  Error  (%) 

0 

0.0572 

0.0572 

0 

1 

0.0574 

0.0573 

0.17452 

9 

0.0714 

0.0715 

0.13986 

11 

-0.0318 

-0.0317 

0.315457 

20 

-0.0162 

-0.0162 

0 

Table  2.8  Electric  field  £’x(xl07V/m)  in  inhomogeneity  InAs/GaAs(lll)  for  different 

mesh  sizes 


^"^-Alesh  size 
x  (nm) 

100 

120 

Relative  Error  (%) 

0 

0.513 

0.515 

0.38835 

1 

0.506 

0.508 

0.3937 

9 

0.110 

0.111 

0.9009 

11 

0.278 

0.281 

1.06762 

20 

-0.171 

-0.171 

0 

Shown  in  Figures  2.5  and  2.6  are  the  hydrostatic  strain  (%x  +yZ7)  along  both  horizontal 
and  diagonal  lines  in  both  InAs/GaAs  (001)  and  InAs/GaAs  (111).  These  results  are 
similar  to  those  for  the  InAs/GaAs  QWR  with  trapezoidal  shape  [Pan05].  In  general,  the 
elastic  strain  fields  inside  the  QWR  are  much  larger  than  those  in  the  substrate  and  the 
difference  in  the  strain  fields  based  on  the  inhomogeneity  and  inclusion  models  is 
apparent,  particularly  within  the  QWR. 
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Hydrostatic  strain  Hydrostatic  strain 


Figure  2.5  Hydrostatic  strains  (yxx  +y//)  in  a  square  QWR  InAs/GaAs  (001)  along  the 

horizontal  (a)  and  diagonal  (b)  lines. 


Figure  2.6  Hydrostatic  strains  (yxx  +y//)  in  a  square  QWR  InAs/GaAs  (111)  along  the 

horizontal  (a)  and  diagonal  (b)  lines. 
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Ex(X107v/m)  Ex(X107v/m) 


Figure  2.7  Ex  (a)  and  E,  (b)  along  the  horizontal  line,  and  Ex  (c)  and  E,  (d)  along  the 
diagonal  line,  in  square  QWR  InAs/GaAs  (111). 


It  has  been  well-known  that  while  there  is  no  induced  electric  field  in  the  InAs/GaAs 
(001)  system,  large  electric  fields  can  be  observed  in  QWs,  QWRs,  and  QDs  (111) 
systems  [Pan02b,Smi86].  Here  in  Figure  2.7  we  show  for  the  first  time  that  the  induced 
electric  fields  along  the  horizontal  and  diagonal  lines  in  InAs/GaAs  (111)  of  a  square 
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QWR  can  be  large  and  that  the  difference  of  the  electric  fields  based  on  the  inclusion  and 
inhomogeneity  models  can  be  significant,  especially  within  the  QWR.  In  other  words, 
electric  fields  in  the  InAs/GaAs  (111)  orientation  should  not  be  neglected,  and  it  should 
be  considered  using  the  correct  inhomogeneity  model  as  the  simple  inclusion  model 
could  be  completely  wrong.  Furthermore,  for  the  induced  electric  field,  its  magnitudes 
both  inside  and  outside  the  QWR  are  comparable,  in  contrast  to  the  corresponding  strain 
field  featured  in  Figures  2.5  and  2.6. 


CT5 

CO 

o 


c n 
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Figure  2.8  Hydrostatic  strains  (yxx  +yzz)  in  hexagon  QWR  InN/AIN  (0001)  along 
horizontal  (a)  and  inclined  (b)  lines. 


Figure  2.8  shows  the  hydrostatic  strain  in  InN/AlN(0001)  along  the  horizontal  and 
inclined  lines.  Similar  to  the  hydrostatic  strain  in  the  InAs/GaAs  system,  we  observed 
that,  outside  the  QWR,  both  the  inclusion  and  inhomogeneity  models  predict  agreeable 
results,  with  their  magnitudes  outside  being  also  much  smaller  than  those  inside  the 
QWRs.  However,  near  the  interface  between  the  QWR  and  substrate,  apparent  difference 
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can  be  observed.  Furthermore,  for  the  hydrostatic  strain  inside  the  QWR,  the  inclusion 
model  predicts  a  high  value  than  the  inhomogeneity  one  (about  11%).  For  the  InN/AIN 
(1000),  the  hydrostatic  strain  distribution  is  similar. 

Figure  2.9  shows  the  electric  field  distribution  along  the  horizontal  line  for  both 
InN/AIN  (0001)  and  (1000).  For  orientation  (0001),  the  only  non-zero  E-field  is  Ex  whilst 
for  (1000),  we  have  only  Ez.  It  is  observed  from  Fig.  2.9  that  the  inclusion  model  predicts 
very  close  results  as  compared  to  those  based  on  the  inhomogeneity  one,  except  for  the 
points  close  to  the  interface  where  the  E-field  experiences  a  shape  change,  resulting  in 
different  values  based  on  different  models. 

Shown  in  Figures  2.10  and  2.11  are  the  E-fields  Ex  and  E,  along  the  inclined  line  in 
both  InN/AIN  (0001)  and  (1000).  It  is  noted  that  along  the  inclined  line,  both  Ex  and  E, 
are  different  from  zero,  as  compared  to  Figure  2.9  where  one  of  the  electric  components 
is  zero.  Similar  to  Figure  2.9,  the  results  predicted  by  the  simple  inclusion  model  are 
mostly  reliable  as  compared  to  those  based  on  the  inhomogeneity  model;  however,  inside 
the  QWR  or  close  to  the  interface  between  the  QWR  and  substrate,  the  inhomogeneity 
should  be  employed. 
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gure  2.9  Ex  in  InN/AlN(1000)  (a)  and  E,  in  InN/AlN(0001)  (b)  along  the  horizontal  line. 
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Figure  2.10  Ex  (a)  and  E,  (b)  in  InN/AIN  (0001)  along  the  inclined  line. 
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Figure  2.1 1  Ex  (a)  and  E,  (b)  in  InN/AlN(1000)  along  the  inclined  line. 


2.5.3  QWRs  of  Different  Polygonal  Shapes 

In  this  example,  we  study  the  induced  E-fields  inside  and  outside  a  polygonal  QWR 
with  different  sides  (n=3,4,5,6,10,  and  infinity),  a  model  used  before  for  the 
corresponding  comer  singularity  study  [Pan06].  The  model  is  shown  in  Figure  2.12  with 
the  E-fields  being  calculated  along  the  horizontal  x-axis.  Only  the  results  from  the 
inhomogeneity  models  InAs  (lll)/GaAs  (111),  InN  (0001)/ AIN  (0001),  and  InN 
(0001)/A1N  (0001)  are  presented. 
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Figure  2.13  shows  that  inside  and  outside  the  QWR,  both  E-field  components  are 
nonzero  along  the  horizontal  jt-axis  and  that  just  as  for  the  second  example,  the 
magnitudes  of  the  E-fields  are  comparable  both  inside  and  outside  the  QWR.  We  also 
observe  that  the  results  from  the  regular  triangle  and  square  QWRs  are  completely 
different  from  the  other  polygonal  QWRs.  This  can  be  seen  more  clearly  from  Table  2.9 
where  the  E-fields  at  the  center  of  the  polygons  are  listed.  It  is  obvious  that  the  signs  of 
Ex  and  Ez,  are  respectively  the  same  for  both  triangle  and  square  QWRs,  which  have 
opposite  signs  as  compared  to  those  in  other  polygons.  Furthermore,  while  the  central 
values  of  the  E-fields  are  very  small  for  the  triangle  QWR,  the  Ez  component  in  square 
QWR  are  much  larger  than  those  in  other  polygon  (50-100%  large). 
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Figure  2.13  Ex  (a)  and  E,  (b)  in  InAs/GaAs  (111)  along  the  horizontal  line  for 
n=3,4,5,6,10,  and  infinity  (circle). 


Figure  2.14  show  the  E-field  along  the  horizontal  axis  in  InN/AIN  (1000)  and  (0001). 
Due  to  the  symmetric  property  of  the  problem  the  only  nonzero  component  is  Ex  in  (1000) 
orientation  and  Ez  in  (0001)  orientation.  It  is  noted  that: 
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(1) .  There  is  a  sharp  change  in  the  E-field  at  the  geometric  corner  point  (x=10nm); 

(2) .  The  trend  of  the  field  variation  in  triangle  QWR  is  completely  different  to  those  in 
other  polygons; 

(3) .  At  the  center  of  the  polygons,  the  E-field  of  the  square  is  much  larger  than  those  in 
other  polygons  where  the  results  are  very  close  to  each  other  (Table  2.10)  except  for  the 
square  QWR.  The  difference  of  the  E-field  magnitude  at  the  center  between  square  and 
other  polygons  is  about  25%. 


Table  2.9  Electric  fields  at  the  center  of  the  polygonal  QWRs  of  InAs/GaAs  (1 1 1) 


Polygon 

sides/Fields 

Ex  (MV/m) 

Ey  (MV/m) 

3 

-0.0071 

0.0012 

4 

-0.1701 

1.1285 

5 

0.2182 

-0.1546 

6 

0.2396 

-0.1430 

10 

0.2613 

-0.1851 

Circle 

0.2689 

-0.1902 

Table  2.10  Electric  fields  (xlOv/m)  at  the  center  of  the  polygonal  QWRs  of  InN/AIN 


Polygon 

sides/Fields 

Ex  in  (1000) 

Ez  in  (0001) 

3 

-42.33 

-45.74 

4 

-56.89 

-56.89 

5 

-42.86 

-43.00 

6 

-41.73 

-43.96 

10 

-42.63 

-42.76 

Circle 

-42.66 

-42.66 
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Figure  2.14  Ex  (a)  in  InN/AIN  (1000)  and  Ez  (b)  in  InN/AIN  (0001)  along  the  horizontal 

line  for  n=3,4,5,6,10,  and  infinity  (circle). 
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2.6  Conclusions 


In  this  chapter,  an  accurate  BEM  modeling  is  proposed  for  the  strain  and  electric  fields 
analysis  in  QWR  structures.  Constant  elements  are  employed  to  discretize  the  interface 
between  the  QWR  and  substrate  with  the  integrand  being  analytically  carried  out  utilizing 
the  exact  close-form  Green’s  functions  for  anisotropic  piezoelectric  solids.  The  material 
models  are  both  group  III-V  and  nitride  III  group,  and  both  the  well-accepted  inclusion 
and  inhomogeneity  models  are  discussed.  The  elastic  and  piezoelectric  response  at  any 
location  can  be  predicted  based  on  the  inclusion  and  inhomogeneity  models.  From  our 
study,  some  important  features  are  observed,  with  those  for  the  strain  field  being 
consistent  to  recent  published  results:  1).  In  the  substrate  and  far  away  from  the  QWR, 
both  the  inclusion  and  inhomogeneity  models  predict  very  close  strain  field.  In  other 
words,  the  simplified  homogeneous  inclusion  model  can  be  safely  applied  if  one  could 
like  to  have  a  quick  estimation  of  the  strain  field.  2).  For  points  inside  or  close  to  the 
QWR,  the  strain  difference  between  the  two  models  can  be  as  high  as  10%  for  the  test 
structures,  which  can  result  in  strong  variations  of  the  confined  electronic  states.  3). 
While  the  magnitude  of  strain  inside  the  QWR  is  much  larger  than  that  outside,  the 
electric  fields  have  the  same  magnitudes  both  inside  and  outside  the  QWR.  4).  While  the 
difference  in  the  electric  fields  based  on  both  inclusion  and  inhomogeneity  models  is 
small  in  the  InN/AIN  system,  that  in  the  InAs/GaAs  (111)  can  be  large.  In  other  words, 
the  inhomogeneity  model  should  be  employed  if  possible.  5).  It  is  also  observed  that  the 
electric  fields  in  the  QWR  depend  strongly  on  the  QWR  geometry  shape;  The  electric 
field  in  triangle  and  square  QWRs  is  different  from  those  in  the  polygons  made  of  more 
than  4  sides. 
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CHAPTER  III 


ELASTIC  AND  ELECTRIC  FIELDS  INDUCED  BY  QWR  IN  BIMATERIAL  PLANE 

3.1  Summary 

Two-dimensional  BEM  is  very  convenient  in  analyzing  stress  and  fracture  problems  in 
bimaterial  solids.  The  bimaterial  Green’s  functions  in  transversely  isotropic  piezoelectric 
solids  [Din04],  anisotropic  elastic  [Pan99a]  and  piezoelectric  bodies  [Pan99b,  Gar05, 
Gro05],  and  in  magnetoelectroelastic  solids  [Qin05]  were  studied  and  applied  to  different 
mechanical  and  piezoelectric  problems.  We  remark  that  these  Green’s  functions  were 
presented  in  the  Lekhnitskii  formalism  [Pan99a].  Applications  of  the  2D  BEM  also 
include  the  determination  of  stress  intensity  factors  for  interfacial  rigid  line  inclusion 
[Lee02],  symmetric  BEM  formulation  for  cohesive  interface  problem  [Sal03],  and 
analysis  of  collinear  interfacial  cracks  interaction  [Sun06].  Other  extended  studies 
include  thermomechanical  closure  of  interfacial  cracks  [Gia07],  fracture  problems  in 
magnetoelectroelastic  composites  [Gar07],  and  dynamic  fracture  problems  involving 
bimaterial  solids  [Lei03,  Lei04,  Lei07]. 

In  this  chapter,  we  develop  the  anisotropic  bimaterial  Green’s  functions  in  terms  of  the 
efficient  and  powerful  Stroh  formalism  [Tin96]  and  the  corresponding  BEM  program  for 
the  analysis  of  the  strain  energy  density  and  the  relative  strain  energy  for  a  QWR  free¬ 
standing  on  or  embedded  in  an  anisotropic  substrate.  The  boundary  of  QWR  was 
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discretized  with  constant  elements  so  that  the  involved  kernel  integration  can  be  carried 
out  in  an  exact  closed  form,  which  is  solved  for  the  boundary  (interface)  values.  The 
induced  elastic  field  inside  and  outside  the  QWR  can  be  obtained  using  the  solved 
interface  values.  We  remark  that  in  order  to  apply  our  bimaterial  program  to  the  free¬ 
standing  QWR  case,  we  only  need  to  assume  that  the  material  stiffness  in  the  upper  half¬ 
plane  is  much  smaller  than  that  in  the  lower  half-plane  substrate. 

Our  bimaterial  BEM  program  is  first  tested  against  various  existing  results.  It  is  then 
applied  to  calculate  the  strain  energy  density  and  the  relative  strain  energy  in  the  InAs 
(1 1 1)  QWR  which  is  free-standing  on  or  embedded  in  the  GaAs  substrate. 

This  chapter  is  organized  as  follows:  In  Section  3.2,  the  basic  problem  is  presented.  In 
Section  3.3,  the  boundary  integral  equation  and  the  required  four  sets  of  Green’s 
functions  for  the  bimaterial  matrix  are  presented.  In  Section  3.4  and  3.5,  the  integration 
and  derivatives  of  bimaterial  Green’s  function  are  obtained.  In  Section  3.6,  the 
formulation  of  elastic  strain  energy  density  and  the  relative  strain  energy  are  presented. 
Numerical  examples  are  presented  in  Section  3.7,  and  conclusions  are  drawn  in  Section 
3.8. 

3.2  Problem  Description  for  QWR  in  Bimaterial  Substrate 

We  assume  that  there  is  a  QWR  of  arbitrary  shape  within  a  bimaterial  substrate.  A 
uniform  eigenstrain  y*j  field  is  applied  inside  the  QWR  as  shown  in  Figure  3.1.  The 

elastic  moduli  in  the  QWR,  material  1  and  material  2  are  denoted,  respectively,  as  c]'jkl , 
cL  and  cL .  In  each  domain,  the  basic  equations  remain  the  same  as  in  section  2.2. 
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Figure  3.1  A  general  QWR  within  a  bimaterial  substrate:  An  eigenstrain  is  applied  to 
the  QWR  which  is  an  arbitrarily  shaped  polygon. 

3.3  Bimaterial  Green’s  Function 

The  boundary  of  QWR  is  discretized  into  constant  elements.  The  discretized  integral 
equations  remain  same  as  Eqs.  (2.20)  and  (2.21).  It  is  noted  that  in  order  to  solve  the 
discretized  boundary  integral  equations,  the  involved  Green’s  functions  f/jj  and  Tjj  and 
their  integrations  on  each  element  are  required.  For  the  QWR  domain,  only  the  full-plane 
Green’s  functions  are  required,  which  can  be  expressed  as  (2.22)  and  (2.23). 

For  the  matrix  domain  (i.e.,  the  bimaterial  substrate),  due  to  the  relative  locations  of 
the  source  and  field  points,  there  are  four  sets  of  Green’s  functions.  Again,  these  Green’s 
functions  were  derived  in  terms  of  the  Lekhnitskii  formalism  [Pan99a].  However,  the 
Stroh  formulism  is  more  convenient  and  efficient  [Tin96]  because  the  it  possesses  the 
eigenrelation  that  relates  the  eigenmodes  of  stress  functions  and  displacement  to  the 
material  properties.  Therefore,  we  first  present  these  bimaterial  Green’s  functions  in  the 
Stroh  formulism.  We  refer  to  Figure  3.1  where  materials  1  and  2  occupy  the  half-plane 
z>0  and  z<0,  respectively. 
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First,  let  the  source  point  X  =(X,  Z)  be  in  material  1  (z>0).  Then,  if  the  field  point  x=(x, 
z)  is  in  the  z>0  half-plane,  the  displacement  and  traction  Green’s  functions  can  be 
expressed  as 


u\  =-Im{4  In  (zr-sr)Al  +  2K  Hz, 

ft  V=1 


K  =-1imi b),  ai +j}B\  p'n‘  zh  e;;-*]) 

ft  Z,.-sr  Zr  ~  s 


On  the  other  hand,  if  the  field  point  is  in  the  z<0  half-plane,  then 


n 


V— 1 


(3.1) 

(3.2) 


(3.3) 


n  ~r  ■ 


2  PrH  1  n3  q!2,v j 


(3.4) 


Now,  let  the  source  point  (X,  Z)  be  in  material  2  (z<0).  Then,  if  the  field  point  (x,  z)  is 
in  the  z>0  half-plane,  we  have 


=iimyjA>u, -.oaf] 


kj  jr 

v=l 


T!l=-Umf}B)rl’'n'  n'Qf  ] 
^  7^  zr-sv 


On  the  other  hand,  if  the  field  point  is  in  the  z<0  half-plane, 

u  l  =  -  Im{  A),  ln(zr  -  ,,  )A;,  +  Hz,-  J,  )&  ]  | 


K 


V— 1 


71  Z~  .V„  Z~  S„ 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


In  Eqs.  (3.1)-(3.8),  the  superscript  1  or  2  to  p,  A,  and  B ,  denotes  that  these  eigenvalues 
and  the  corresponding  eigenmatrices  belong  to  material  1  or  2.  The  matrix  Q  depends  on 
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the  relative  locations  of  the  source  and  field  points.  Assuming  that  the  source  point  is  in 
the  half-plane  of  material  2  (1=1  or  2).  Then  Q  can  be  expressed  as 

Qu,v  =  Ku{iv)jjA  (3.9) 

if  the  field  point  is  in  the  half-plane  of  material  2,  and 

Qlk"  =  Krfd,),K  (3-10) 

if  the  field  point  is  in  the  other  half-plane  of  material  //  (ju  ±  2).  In  Eqs.  (3.9)  and  (3.10), 
the  matrix  K  is  given  by 

Ku  =(AAyl(MA  -Ma)Aa  (3.11) 

kAm  =  c AMy\MM  +ma)~\ma  +  ma)aa  (3.12) 

with  Ma  =-iBA(AAy1  (  2=1  or  2)  and  the  diagonal  matrix  Iv  has  the  following 
expression  for  different  indexes  v, 

/t=diag[  1,0,0];  /2=diag [0,1,0];  /3==diag[0,0,l]  (3.13) 

3.4  Integration  of  Bimaterial  Green’s  Function 

In  order  to  carry  out  the  line  integral  of  the  involved  Green’s  functions  over  each 
constant  element  (e.g.,  along  the  n-th  element  Th),  we  let  the  generic  element  start  from 

point  1  (xi,  zi)  and  end  at  point  2  (x2,  z2),  with  length  /  =  yj(x2  -  x^2  +{z2  ~  z^2  ■  Then, 
the  integration  of  the  involved  Green’s  displacements  and  tractions  (only  the  integral 
variable-dependent  function)  for  the  QWR  (which  requires  only  the  infinite-plane 
Green’s  functions)  can  be  carried  out,  which  has  been  expressed  in  section  2.4. 
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We  point  out  that  similar  expressions  can  be  found  for  the  integration  of  the  bimaterial 


Green’s  functions  in  the  substrate.  Let  the  source  point  X  be  in  material  1,  if  the  field 
point  x  is  in  material  1,  the  corresponding  integrals  are 

K(x,z)  =  \r  In (z\  ~s\)dTn 


=  /[ 


("  +  /;-)-  S  '  inl^4^  V'  1  +  ln(x2  +  p\z2  -  - 1] 


(3.14) 


w 


1  AAAL  j 

(x2  -X,)  +  pr  (z2  -  Z1 )  (x,  +  PrZ\  )  -  s 

'L  (x,  z)  =  Jr  i n (> '  -5,‘y/r, 


=  /[  (xi  + pUi)~s! 


ln[ 


{x2  + p].z2)-sl 


- i - AAA L - 1 - ~i 

(x2-x1)  +  pr(z2-z1)  ( xl+prz1)-sv 


]  +  ln(x2  +  p\z2- s])  - 1] 


(3.15) 


and 


g 


n  Zl  —  S  1 


-dr„=i[ 


ln(^4^4)]  (3-16) 


(x2  -XJ  +  p1r(z2-z1)  x,  +  p1rz1  -s 


dlrv(x,z)  =  f  ,  1  _tdTn  =  /[ 

•'in  -7  —  e 


ln(*2+/V2  4)J  (3-17) 


(x2  -xl)  +  pl(z2  -  Z[)  Xy  +  p^-s] 


When  source  point  X  is  in  material  1,  and  the  field  point  x  is  in  material  2,  the 
corresponding  integrals  are 


hr!  (*.  z)  =  Jr  In (z;  -  si  )dTn 


/[ 


(Xj  +P2rzl)-s\ 


Xx2+  p;z2)~  si 


2 ,  ,  M72  2;  n + in(x2  +  P;z2  -  si )  - 1] 


(3.18) 


(x2-x1)  +  pr(z2  -Z.t)  (Xj  +prZ1)-Sv 


and 


8r 


21 


(^,z)=Jr 


r^n  =  /[ 


1 


ln( 


X2+PrZ2  ~Sv 


(x2-xl)  +  pr(z2-zl)  Xx  +  p;.z 


)]  (3.19) 


When  source  point  X  is  in  material  2,  and  the  field  point  x  is  in  material  1,  the 
corresponding  integrals  are 
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(3.20) 


h); (x, z)  =  jr  Inf-)  -,s 
(*1  +P\Z1)~SI 


wn 


=  /[ 


ln[ 


(x2  +  p].z2)  -  s; 


(x2  -  X, )  +  pi  (z2  -  Zl )  (x,  +  p].zl )  -  5 


J  +  ln(x2  +  p\z,2-s2v)-\\ 


and 


8^  (x,  z)  =  f  ,  1  2drn  =  /[ 
Jr»  z,.  —  s„ 


1 


\n(X2+Pf2  %]  (3-21) 


(x2  -X^  +  p^Zj-Zj)  Xj  +/7rz,  -5 


When  source  point  X  is  in  material  2,  and  the  field  point  x  is  in  material  2,  the 
corresponding  integrals  are 

h; (x,  z)  =  Jr  ln(-;  -  s; )dTn 


and 


/[ 


(x,  +P2rzx)-s 


ln[ 


(x2  +  p;z2)-s; 


]  +  ln(x2  +  p1rz1-s2r)-\\ 


/  \  .  2  /  \L/  .  2  \  2  J  \  2  *  r  2 

(x2  -  Xj )  +  P,.  (z2  -  Z, )  (Xj  +  Pr  Zl)-Sr 


w. 


(x,z)  =  l  \n(z2-s2)dYn 


=  /[ 


(x,  +p;z1)-5l2 


(x2  -Xj)  +  pr(z2  -Zj)  (Xj+^zJ-s 


ln|IJj  +  ;’rjl  +  Infc  +  )  - 1] 


1 


X2  +  p;z2 


(3.22) 


(3.23) 


- 2 - Inf— - - - 

(x2-x1)  +  pr(z2-z1)  Xj  +  pr  Zj  —  5 


-)]  (3.24) 


d. 


-  z  .  z 

-z;-sv  (x2-x1)  +  pr(z2-zl)  Xl+prZl~S 


2/  +  (3.25) 


3.5  Derivatives  of  Bimaterial  Green’s  Function 

In  order  to  obtain  the  strain  and  stress  fields,  we  also  need  to  take  the  derivative  of  the 
Green’s  displacement  and  traction  with  respect  to  the  source  point,  and  then  find  the 
corresponding  integration.  The  integration  of  the  involved  Green’s  function  derivatives 
for  the  QWR  domain  can  be  written  as 
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V."Jr 


3ln(z,,-5f) 

c)x 


clT„  =  — 


/ 


ln( 


A  +  prz2  -  sr 


(x2  -x,)  +  pr  (z2  -  Z,  )  X,+  PrZI  -  s 


)  (3.26) 


h  =  p  h 

r,z  *  r  r,x 


(3.27) 


*~-L 


[ - - - + - - - ]  (3.28) 


dx  zr  —  sr  "  (x2-x1)  +  pr(z2-zl)  x2  +  prz2-sr  Xx  +  przx-s 


Sr.z  Pr8r,x 


(3.29) 


Again,  similar  expressions  can  be  found  for  the  corresponding  bimaterial  substrate. 


3.6  Strain  Energy  Density 

During  growth  of  QWR  nanostructures,  strain  energy  in  the  QWR  plays  an  important 
role.  Therefore,  its  calculation  and  prediction  are  of  great  interest.  Following  [Mur87], 
the  strain  energy  in  the  QWR  can  be  expressed  as 


We  remark  that  the  integrand 


w=jlr<7iJ(rt-r;>dv 


(3.30) 


is  the  strain  energy  density  within  the 


QWR,  which  will  be  also  calculated  and  studied. 

Applying  the  Gauss  divergence  theorem  and  assuming  that  the  eigenstrain  is  uniform 
inside  the  QWR,  the  total  elastic  strain  energy  of  QWR  in  Eq.  (3.30)  can  be  expressed 
alternatively  as 


(3.31) 


where  T  is  the  boundary  of  QWR.  While  the  first  boundary  integration  on  the  right-hand 
side  of  Eq.  (3.32)  represents  the  work  done  by  the  traction  acting  on  the  boundary,  the 
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second  boundary  integration  represents  the  work  done  by  the  traction  associated  with  the 
eigenstrain  as  we  recall  that  t^'  =Cijkly*jn] .  The  last  volumetric  integration  corresponds 

to  the  initial  constant  strain  energy  of  the  eigenstrain  field,  i.e.,  the  constant  strain  energy 
in  the  wetting  layer  (i.e.,  the  thin  layer  covers  the  substrate)  due  to  the  misfit  strain 
[Hol99], 

In  QWR  growth,  one  of  the  important  parameters  is  the  relative  strain  energy  A 
[Hol99,  Sta04].  Namely,  the  ratio  of  the  strain  energy  change  relative  to  the  constant 
strain  energy  in  the  wetting  layer  due  to  the  misfit  strain  over  the  constant  strain  energy, 
as  defined  below 

W -  if  c^rldv  \UutdT-\\  t[>-\dr 
A  =  2  - =  -2—— - -  (3.32) 

cmr,fddv  -j Vcliur;fudv 

where  A  is  also  called  the  relative  strain  energy,  which  will  be  numerically  examined  in 
the  next  section.  We  remark  that  while  various  energetic  parameters  were  introduced  for 
characterizing  defect  dynamics  [PhiOl],  this  relative  strain  energy  has  been  successfully 
applied  in  the  prediction  of  new  QWR  formation  and  QWR  array  patterns  under  the 
Stranski-Krastanow  growth  mode  [Hol99,  Sta04,  Yan03b]. 

3.7  Numerical  Examples 

Before  applying  our  solution  to  examine  the  strain  energy  in  QWR  system,  we  have 
first  checked  our  program  for  the  reduced  cases  with  existing  solutions  [Pan04b].  The  test 
example  is  InAs  QWR  with  size  20nm  by  20nm  embedded  in  half-space  GaAs  substrate. 
The  boundary  condition  is  assumed  to  be  traction  free  on  the  half-space  surface.  The 
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strain  field  along  the  dash  line  in  Fig.  3.2  are  calculated  and  compared  with  results  in 
[Pan04b].  From  Fig.  3.3,  we  can  see  the  two  solutions  match  each  other  very  well.  Other 
numerical  tests  have  also  carried  out,  all  showing  that  our  solution  and  program  are 
correct. 


z 


Figure  3.2  A  square  InAs  QWR  embedded  in  GaAs  substrate. 
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Figure  3.3  Strain  comparison  of  bimaterial  solution  and  analytical  solution  (a)yxx  (b)yzz. 
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Therefore,  after  validation,  we  now  apply  our  solution  to  calculate  the  relative  strain 
energy  and  the  distribution  of  the  strain  energy  density.  We  remark  that  while  our  BEM 
formulation  can  be  applied  to  more  complicated  situations,  here  we  only  study  the 
reduced  case  where  the  QWR  is  free-standing  on  or  embedded  in  the  half-plane  substrate. 
The  half-plane  model  is  reduced  from  our  bimaterial  result  by  setting  the  upper  half-plane 
(z>0)  with  very  low  material  stiffness  as  compared  to  the  lower  half-plane  substrate  (The 
material  stiffness  in  the  upper  half-space  is  equal  to  10"10  times  the  material  stiffness  in 
the  lower  half-space). 

3.7.1  Variation  of  Relative  Strain  Energy  with  QWR  Depth 

We  first  apply  our  bimaterial  BEM  program  to  study  the  problem  where  a  square 
InAs(lll)  QWR  growing  on  a  GaAs(lll)  substrate  (Figure  3.4).  The  material  properties 
of  InAs  and  GaAs  in  the  (11  Indirection  are  obtained  by  coordinate  transform  from  those 
in  the  (001) -directions  [Pan07].  The  QWR  has  a  dimension  of  20nmx20nm  and  a  uniform 
misfit  strain  field  y'xx  =  y*  =  y*^  =0.07.  The  boundary  condition  on  the  surface  of  the 

substrate  is  assumed  to  be  traction-free.  While  a  total  free-standing  QWR  on  the  substrate 
is  illustrated  in  Figure  3.4a,  a  fully  embedded  QWR  is  shown  in  Figure  3.4c.  We  let  cl  be 
the  depth  of  the  QWR,  measured  from  the  bottom  side  of  the  QWR  to  the  surface  of 
substrate  (Figure  3.4b).  We  then  calculate  the  relative  strain  energy  (i.e.,  Eq.  3.32)  of  the 
QWR  as  a  function  of  depth  d,  varying  from  -lOnm  (Figure  3.4a)  to  70nm,  at  the  interval 
of  5nm.  Figure  3.5  shows  clearly  that  with  increasing  depth  d,  the  magnitude  of  the 
relative  strain  energy  increases,  reaching  the  value  when  the  QWR  is  within  an  infinite 
substrate  (=-0.19).  In  other  words,  bringing  an  embedded  QWR  to  the  surface  will 
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decrease  the  relative  strain  energy.  This  is  true  since  with  increasing  surface  area  of  the 
QWR  to  the  air,  more  surface  energy  will  be  released,  resulting  in  small  relative  strain 
energy.  This  important  feature  could  represent  the  competition  between  the  surface  and 
bulk  energies  i.e.  [Pao06a,  Pao06b],  and  should  be  particularly  interesting  to  epitaxial 
growth. 


Figure  3.4  Geometry  of  a  square  InAs  QWR  on  or  inside  the  GaAs  substrate.  Shown  in  (a) 
to  (c)  are  the  three  special  cases:  Namely,  totally  free-standing  (a),  half  free-standing 
(half-in  and  half-out  in  (b)),  and  fully  embedded  (c). 


Figure  3.5  Variation  of  relative  strain  energy  with  depth  of  a  square  QWR. 


51 


3.7.2  Effect  of  QWR  Shape  on  Relative  Strain  Energy  and  Distribution  of  Strain  Energy 
Density 

We  assume  now  that  there  is  an  isosceles  triangle  of  InAs  (111)  QWR,  free-standing 
on  the  GaAs  (111)  substrate.  The  base  angle  of  the  triangle  varies  from  30°  to  75°  while 
the  area  of  triangle  maintains  the  same.  The  misfit  strain  is  again  uniform,  i.e.  y*xx  = 

y*v  =  y  =0.07.  The  boundary  condition  on  the  surface  of  the  substrate  is  traction-free. 

Both  the  relative  strain  energy  and  strain  energy  density  in  the  QWR  are  investigated. 

The  relative  strain  energy  for  the  triangular  QWR  with  different  base  angles  are  listed 

in  Table  3.1  and  its  variation  with  the  base  angle  are  also  shown  in  Figure  3.6.  It  is 

observed  that  the  magnitude  of  the  relative  strain  energy  increases  with  increasing  base 

angle  of  the  QWR  triangle.  In  other  words,  the  steeper  the  QWR  is,  the  larger  the 

magnitude  of  the  relative  strain  energy  becomes. 

Table  3.1  Variation  of  relative  strain  energy  with  base  angle  of  the  isosceles  triangle  of 
InAs  (111)  QWR,  which  is  free-standing  on  the  GaAs  (111)  substrate. 


Base  angle 
(degree) 

Relative  strain  energy 

30 

-5.88E-03 

37 

-6.91E-03 

45 

-7.90E-03 

52 

-8.70E-03 

60 

-9.76E-03 

67 

-1.11E-02 

75 

-1.43E-02 

52 


Figure  3.6  Variation  of  relative  strain  energy  with  base  angle  of  isosceles  triangle  of  InAs 
(111)  QWR,  which  is  free-standing  on  the  GaAs  (111)  substrate. 

The  contours  of  the  strain  energy  density  for  these  free-standing  triangular  QWRs  are 
plotted  in  Figure  3.7.  It  is  clear  that  the  strain  energy  density  distribution  is  strongly 
influenced  by  the  QWR  shape.  Particularly  at  the  locations  near  the  vertex  and  base 
corners,  the  strain  energy  densities  are  quite  different  among  these  triangles.  We  further 
observed  that  the  strain  energy  density  near  the  vertex  is  larger  than  those  near  the  base 
corners.  The  magnitude  of  strain  energy  density  increases  roughly  from  2.15  to  2.35 
(109kg/ms2)  as  the  base  angle  varies  from  30°  to  75°. 


53 


2.15 

2.1 

2.05 

2 

1 .95 
1 .9 
1 .85 
1 .8 
1 .75 
1  .7 
1 .65 
1 .6 
1  .55 
1 .5 
1  .45 
1 .4 
1  .35 
1  .3 
1 .25 
1  .2 


1 


2.2 

2.1 

2 

1.9 

1.8 

1.7 

1.6 

1.5 

1.4 

1.3 

1.2 

1.1 


I 


2.3 
2.2 
2.1 
2 

1.9 

1.8 

1.7 

1.6 

1.5 

1.4 
1.3 
1.2 
1.1 


54 


Normalized  x 
;y  density  (10 

(111)  QWR  with  base  angle  30°(a),  45°(b),  60°(c),  and  75°(d). 


9  2 

Figure  3.7  Contours  of  strain  energy  density  (10  kg/ms  )  in  isosceles  triangle  of  InAs 


3.8  Conclusions 

In  the  chapter,  we  first  derive  the  bimaterial  Green’s  functions  in  anisotropic  elastic 
media  in  terms  of  the  elegant  Stroh  formalism.  The  corresponding  BEM  formulation  is 
then  presented.  Since  the  involved  Green’s  functions  are  in  exact  closed  forms,  the  kernel 
integration  can  be  analytically  carried  out  for  the  constant  element  discretization.  After 
testing  our  bimaterial  BEM  program  for  various  reduced  simple  cases,  we  then  apply  our 
program  to  calculate  the  strain  energy  density  and  the  relative  strain  energy  in  InAs  (111) 
QWR  free-standing  on  or  embedded  in  GaAs  (111)  substrate.  Our  numerical  results 
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showed  that,  for  the  case  of  an  isosceles  triangle  of  InAs  (111)  QWR  on  the  substrate,  the 
strain  energy  density  within  the  QWR  strongly  depends  on  the  base  angle  of  the  triangle. 
We  also  observed  that  the  magnitude  of  the  relative  strain  energy  for  this  case  increases 
with  increasing  base  angle.  For  a  square  QWR  either  free  standing  on  or  embedded  in  the 
substrate,  we  found  that  the  magnitude  of  the  relative  strain  energy  increases  with 
increasing  depth  of  the  QWR  in  the  substrate,  which  represents  the  competition  between 
the  surface  and  bulk  energies.  Our  numerical  examples  also  demonstrate  the  significant 
influence  of  the  QWR  shape  and  location  on  both  the  strain  energy  density  and  the 
relative  strain  energy  in  the  QWR,  which  should  be  particularly  useful  to  the  successful 
growth  of  QWRs  via  epitaxial  approach. 


56 


CHAPTER  IV 


STAIN  ENERGY  ON  THE  SURFACE  OF  HALF-SPACE  SUBSTRATE  INDUCED 

BY  AN  ISOLATED  QD 


4.1  Summary 

Novel  superlattice-based  quantum  devices  have  attracted  considerable  attention  in 
recent  years  [Bim98,  Sta04,  Jal05;  Fer07]  where  quantum  dots,  nanowires  and  nanorods 
have  been  intensively  investigated  [Leo93,  Gos95,  Bax05].  Self-assembled  quantum  dots 
have  been  of  particular  interest  because  of  their  three-dimensional  confinement  of  charge 
carriers  and  excitons  and  hence  their  potential  in  nano-scale  electronic  devices  [Pet03, 
Fri07,  Lev07]. 

It  is  well  known  that  the  elastic  fields  produced  by  QDs  substantially  influence  surface 
feature  formations  that  consequently  affect  the  electronic  band  structure.  As  either  a 
detector  or  emitter  of  electromagnetic  signals,  the  bandgap  energy  determines  the 
operational  wavelength  in  say,  the  photovoltaic  mode,  of  the  element.  Hence,  the  elastic 
fields  induced  by  QDs  have  to  be  studied  in  order  to  obtain  a  well  ordered  QD  structure. 
To  analyze  the  elastic  fields  in  and  around  the  QDs  and  other  quantum  structures,  several 
methods  have  been  proposed,  such  as  the  finite  element  and  finite  difference  methods, 
atomistic  modeling,  and  analytical  approaches.  In  particular,  the  Green’s  function 
solutions  have  been  proposed  and  applied  to  QD  studies  [PanOl)].  Because  of  their 
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accuracy  and  efficiency,  the  analytical  methods,  particularly  the  Green’s  function  method, 
are  arguably  more  appealing  to  experimentalists  and  physical  device  designers  in  the 
studies  of  QD  structures  [And99,  PeaOO].  More  recently,  the  multiscale  Green’s  functions 
have  been  also  introduced  to  model  QD-induced  strain  field  [Tew04a,  Tew04b,  Yan05, 
Rea07],  which  have  been  shown  to  be  computationally  efficient.  However,  while  the 
lattice  Green’s  functions  have  been  well-developed  so  far,  the  corresponding  continuum 
Green’s  functions  still  require  further  study,  particularly  when  considering  semiconductor 
anisotropy. 

In  this  chapter,  we  present  an  analytical  method  for  the  QD-induced  strain  field  in 
half-space  semiconductor  substrates  under  the  assumption  of  continuum  elasticity. 
Furthermore,  under  the  epitaxial  growth,  the  misfit  strain  within  the  QD  could  be  gradient, 
instead  of  uniform  distribution  [Mal03],  and  as  such,  the  material  property  in  the  QD 
could  be  difficult  to  calibrate  [Zhu06].  Thus,  the  well-know  inclusion  model  will  be 
adopted  in  this  chapter.  It  was  verified  recently  that  the  inclusion  model  could  predict 
slightly  different  results  as  compared  to  the  inhomogeneity  model  using  the  bulk  property 
of  the  QD  (about  10%  within  the  strained  quantum  structures,  see,  e.g.,  [Yan02,  Pan05]). 
Under  these  assumptions,  we  derive  our  solution  based  on  the  Green’s  function  method  in 
terms  of  the  Stroh  formalism  with  the  corresponding  exact  integration  of  the  Green’s 
functions  over  the  QD  surface  (composed  of  piece-wise  flat  triangles). 

This  chapter  is  organized  as  follows:  In  section  4.2,  the  strained  QD  system  is 
described.  In  section  4.3,  the  surface  of  the  QD  is  approximated  by  a  number  of  flat 
triangles  over  which  the  area  integration  is  carried  out  exactly  so  that  the  induced  elastic 
fields  can  be  expressed  in  terms  of  a  simple  line  integral  over  [0,  Jt[.  If  the  QD  is  a  point 
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source,  then  the  QD-induced  elastic  fields  can  be  analytically  expressed  by  point-force 
Green’s  functions  in  the  half-space.  In  section  4.4,  numerical  examples  are  carried  out  for 
buried  cubic,  pyramidal,  truncated  pyramidal,  and  point  QDs  in  half-spaces.  The  effects 
of  QD  shape  and  depth  on  the  strain  energy  are  discussed.  Conclusions  are  drawn  in 
section  4.5. 

4.2  Problem  description  for  QD  Embedded  in  Anisotropic  Half-space  Substrate 

We  now  assume  that  the  3D  point-force  Green’s  functions  are  given  for  the  half-space 
substrate,  then  for  the  general  eigenstrain  yl  at  x=(xi  ,X2,x?)  within  the  QD  domain  V,  the 

induced  displacement  at  d=(di,d2,di)  can  be  found  using  the  superposition  method.  That 
is,  the  response  is  an  integral,  over  V,  of  the  equivalent  body  force  defined  by  (2.17), 
multiplied  by  the  point-force  Green’s  functions,  as 

ul(d)  =  -fu‘(x;d)[C,l,,,r;jx)lJdV(x)  (4.1) 

where  Uk(x\d )  is  the  j'-th  Green’s  elastic  displacement  at  x  due  to  a  point  force/  point 
charge  in  the  &-th  direction  applied  at  d. 

Integrating  by  parts  and  noting  that  the  eigenstrain  is  nonzero  only  in  the  QD  domain 
V,  Eq.  (4.1)  can  be  expressed  alternatively  as 

uk(d)  =  \ukjXj(x-,d)Cijlmy*m(x)dV(x)  (4.2) 

V 

If  we  further  assume  that  the  eigenstrain  is  constant  within  the  QD  domain  V,  then  the 
domain  integration  can  be  transformed  to  the  boundary  integration.  That  is 

uk (d )  =  Cijlm yl  { U)  (x;d)n,  (x)dS(x)  (4.3) 

ay 
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where  n\(x)  is  the  outward  normal  on  the  boundary  c)V  of  the  QD. 

In  order  to  find  the  elastic  strain,  we  take  the  derivatives  of  Eq.  (4.3)  with  respect  to 
the  observation  point  d  (i.e.,  the  source  point  of  the  point-force/point  charge  Green’s 
function),  which  yields  (k,p=  1,2,3) 

K  a )  =  (  rlc„m  | (*;</)  +  U’Jt  (x ;«/)]*,  (x)dS(x)  (4.4) 

1  sv 

The  stresses  inside  and  outside  the  QD  are  obtained  from  Eq.  (2.14). 

It  is  obvious  that  in  order  to  solve  the  QD-induced  elastic  field,  the  key  is  to  carry  out 
the  surface  integration  involved  in  (4.3)  and  (4.4).  This  requires  the  integral  of  the 
corresponding  half-space  Green’s  functions. 

4.3  Integration  of  Half-space  Green’s  Functions  over  QD  Surfaces 

We  assume  that  the  boundary  of  the  QD  can  be  approximated  by  a  number  of  flat 
triangles.  We  want  to  analytically  integrate  the  half-space  Green’s  functions  over  one  of 
the  flat  triangles.  To  do  so,  we  first  briefly  review  the  half-space  Green’s  functions  in 
general  anisotropic  semiconductors. 

The  half-space  point-force  Green’s  function  with  source  point  at  d  and  field  point  at  x 
can  be  expressed  as  a  sum  of  an  explicit  infinite-space  solution  and  a  complementary  part 
in  terms  of  a  line  integral  over  [0,7t]  [Pan  02b] 

U(x;d)  =  Uoa(x;d)  +  -^-TrAGlATd0  (4.5) 

In2  Jo 

where  overbar  means  complex  conjugate,  superscript  T  denotes  matrix  transpose,  and 
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(4.6) 


(G,),=^r 


(BB)* 


-  p,x3  +  p  d3  - [(x,  -d{) cos G  +  (x2  -d2) sin G ] 


We  point  out  that,  on  the  right-hand  side  of  Eq.  (4.5),  the  first  term  corresponds  to  the 
Green’s  displacement  tensor  in  an  anisotropic  infinite  space.  Its  integration  over  a  flat 
triangle  was  presented  by  Wang  et  al.  [Wan06].  Also  in  Eqs.  (4.5)  and  (4.6),  the  Stroll 
eigenvalues  p3,  and  eigenmatrices  A  and  B  are  all  functions  of  G \  as  well  as  the  elastic 
stiffness  tensor  of  the  materials.  p3  and  A  are  the  solutions  of  the  following  eigenrelation 

[Q  +  p(R  +  RT)  +  p2T]a  =  0  (4.7) 

which  is  exact  same  as  Eq.  (2.24),  but  with 


Qik  -  Ciakpnanp  ’  Rik  -  Ciak3na ,  Tik  -  C,3 


where 


(n.,n2)  =  (cos  G,  sin  G) 


and  a  and  [i  taking  the  values  of  1  and  2.  We  noted  that,  due  to  the  positive  requirement 
on  the  strain  energy  density,  the  eigenvalues  of  Eq.  (4.7)  are  either  complex  or  purely 
imaginary. 

In  order  to  find  the  misfit  strain-induced  elastic  fields,  we  also  need  the  derivatives  of 
the  Green’s  displacement  tensor  with  respect  to  the  field  point  (jq  ,  x2,  x3)  of  the  Green’s 


function.  They  are  found  to  be  (j=  1,  2,  3) 


dU(x-d)  dU°°(x-d )  1  /  \  3 


(4.10) 


where 
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(4.11) 


(G2), 


(*  *)u 

{-PjX 3  +  Pjd3  - [(jc,  —  dy) cos 0  +  (jc2  -r/2)sin#]}2 


(g ,)  =  fifmg[cos#,cos<9,cos#l 

(g2)  =  <^/<3^[sin  <9,  sin  6>,  sin  6^]  (4.12) 

(g3)  =  diag[pl,  p2,  p3] 


Again,  since  the  integration  on  a  flat  triangle  has  already  been  presented  for  the 
infinite  part  of  the  half-space  Green’s  function  [Wan06],  we  only  need  to  present  the 
integration  of  the  complementary  part  of  the  Green’s  function,  i.e.,  the  integration  of  the 


second  term  on  the  right-hand  side  of  Eq.  (4.10). 


Figure  4.1  Sphere  surface  approximated  by  a  number  of  flat  triangles. 

We  consider  first  the  integration  of  the  Green’s  displacement  tensor.  Again,  we 
assume  that  the  QD  surface  can  be  effectively  approximated  by  a  number  of  flat  triangles 
(Fig.  4.1).  Therefore,  the  integral  expression  Eq.  (4.3)  over  a  flat  triangle,  A,  is 

uk  id)  =  Cijlm y]mnl  J U j  (x-d)dA(x)  (4. 13) 

A 
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where  n  is  the  outward  normal  to  A.  Substituting  Eq.  (4.5)  into  (4.13)  and  changing  the 
integration  orders,  the  contribution  from  the  complementary  part  can  be  expressed  as 


uk(d) 


(4.14) 


While  the  outside  line  integration  can  be  easily  carried  out  by  employing  Gaussian 
quadrature,  we  discuss  the  area  integration  over  the  flat  triangle,  which  can  be  done 
analytically  as  will  be  shown  below.  Actually,  since  in  the  expression  for  Gi  in  Eq.  (4.6), 
its  numerator  is  a  function  of  6  only,  the  integration  over  the  flat  triangle  A  needs  to  be 
carried  out  for  the  following  expression  only 


dA(jc) 

-  ptx3  +  pjd3-[(xl  -  dx )  cos  0  +  ( x2  -  d2 )  sin  0 ] 


(4.15) 


Similarly,  in  order  to  find  the  QD-induced  strain  field,  (Eqs.  (4.4)  and  (4.10)),  one 
needs  only  to  carry  out  the  following  area  integration  over  the  flat  triangle. 


dA(jc) 

{ -p,x3  +  pjd 3  -  [(Xj  -  dx ) cos  6  +  (x,  -  d2 )  sin  #]  }2 


(4.16) 


The  area  integration  over  a  flat  triangle  in  Eqs.  (4.15)  and  (4.16)  can  be  carried  out 
exactly.  To  do  so,  we  introduce  the  following  transformation  between  the  global 
coordinate  system  x  (x\,  x2,  *3)  and  local  coordinate  system  £(£i,  £3)  associated  with 

the  flat  triangle  (Figure  4.2) 

*1  ~  -%  a\\  a\2  ai3 

X2~X02  ~  a2l  a22  a23  *22  (4. 17) 

.A  i  .Vq3  r/31  ^3")  6E33  ^^3 


Then,  the  integration  becomes  (n=l  ,2) 
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1 


(4.18) 


where 


h  I h 

Fn(d],d)  =  \d^2  {< 


0  -/,+/,f2/A  Lil'.^7 


[il{dJ,0)^+f2(<d1,0^2+Udj,0)r 


fx{d j,0)  =  -{pta3X  +  an  cos 0  +  a2l  sin 6) 

f2(d j,6>)  =  -( PjCt32  +  an  cos 0+  a22  sin 0) 

f3(d . , 0)  =  —p. (x03  +  a33%3 )  +  Pjd3  +  d{  cos  6  +  d2  sin  0 
-  (x01  +  al3£3)cos6>  -  (x02  +  a23%3) sin  6 


The  integration  can  now  be  carried  out,  and  the  results  are 


FtfpO) 


1  r  Hi  +  A  ln(f2h  +  A  s  ~  Ah  +  A  w  Ah  +  A  xi 
A  fi  -  Ah  / h  Ak  +  A  A  +  AA'h  -fA  +  f 3 


F2(dj,6) 


A 


1 


-ln( 


f2h  +  /3 


)- 


1 


-ln( 


fih  +  A 


A  A  +  fA/h  -AA  +  A  f2~ f A1  h  AA  + A 


)] 


(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 


With  these  exact  expressions,  the  QD-induced  displacement  and  strain  fields  can  be 
finally  expressed  in  terms  of  the  line  integration  over  [0,7t]  (e.g.,  Eq.  (4.14)  for  the 
induced  displacements).  Again  the  line  integration  can  be  carried  out  numerically  using 
8-point  Gaussian  quadrature. 
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Figure  4.2  Geometry  of  the  flat  triangle  A  (with  corners  1,2,3),  and  transformation 
from  the  global  (xi,X2,X3)  to  local  (£i ,^2,^3)  coordinates  where  £3  is  along  the  outward 

normal  direction  of  the  flat  triangle. 

4.4  Numerical  Examples 

We  now  apply  our  analytical  solutions  to  calculate  the  strain  energy  induced  by  a 
buried  QD  within  the  GaAs  half-space  substrate.  The  surface  of  the  substrate  is  traction 
free.  The  QD  is  located  at  depth  d  below  the  surface  (Figure  4.3)  and  the  misfit  strain  is 
hydrostatic,  i.e.  y'xx  =  y*yy  =  y  =  0.07  .  We  point  out  that  this  misfit  strain  is  relatively 

large  and  the  corresponding  nonlinear  influence  will  be  pursued  using  other  methods, 
such  as  the  multiscale  meshless  method  [She04].  As  for  the  QD  shape,  we  assume  it  to  be 
either  cubic,  pyramidal,  truncated  pyramidal,  or  point  type.  The  QDs  have  the  same 
height  /?  (=4nm,  except  for  point  QD).  To  make  all  the  QDs  (including  the  point  QD) 
have  the  same  volume,  we  have  the  base  length  2.155/?  for  cubic  QD,  upper  length  1.79/? 
and  lower  length  2.5/?  for  truncated  pyramid  QD,  base  length  3.732/?  for  pyramid  QD 
(Figure  4.3).  For  the  point  QD  case,  it  is  located  at  the  middle  height  of  the  cubic  QD  (i.e., 


65 


its  vertical  distance  to  the  surface  is  d+h/2).  We  studied  the  effect  of  the  QD  shape  and 


depth  on  the  strain  energy  on  the  surface. 

Surface 


Figure  4.3  Geometry  for  (a)  a  cubic  QD,  (b)  a  truncated  pyramid  QD,  (c)  a  pyramid 
QD,  and  (d)  a  point  QD.  Top  row  is  the  3D  view  and  bottom  row  is  the  vertical  x-z 
plan  view.  All  these  QDs  have  the  same  volume. 

Shown  in  Figures  4.4a-4.4d  are,  respectively,  contours  of  the  normalized  strain  energy 
on  the  surface  of  GaAs  (001)  (top  row)  and  GaAs  (111)  (bottom  row),  induced  by  a 
buried  cubic,  truncated  pyramidal,  pyramidal  or  point  type  QD.  In  this  example,  the  depth 
d= 2nm,  height  h= 4nm,  and  the  strain  energy  is  normalized  by  118.8xl015Nm  (this 
normalization  factor  is  also  used  below  for  the  strain  energy)  (Table  4.1).  It  is  clear  that 
different  QD  shapes  (including  point  QD)  induce  different  strain  energy  distributions  on 
the  surface  of  the  substrate.  Besides  the  difference  on  the  contour  shape,  the  strain  energy 
values  are  also  different.  For  example,  the  contours  with  value  0.3,  0.6  and  0.9 
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(corresponding,  respectively,  to  the  heavy  blue,  red  and  black  curves)  move  towards  the 
center  when  the  QD  shape  changes  from  left  to  right  (i.e.,  cubic,  truncated  pyramid, 
pyramid,  and  point  type).  Table  4.1  lists  the  maximum  strain  energy  values 
corresponding  to  different  QD  shapes  at  different  depths  within  the  GaAs  substrate  with 
both  (001)  and  (111)  orientations.  It  is  clear,  from  Table  4.1,  that  while  the  induced 
maximum  strain  energy  value  decreases  with  increasing  depth,  its  value  on  the  substrate 
GaAs  (001)  is  always  larger  than  that  of  the  corresponding  inclined  substrate  GaAs  (111). 
Furthermore,  the  effect  of  the  QD  shapes  (cubic,  truncated  pyramid,  pyramid,  and  point 
QD)  on  the  strain  energy  is  complicated,  as  can  be  seen  from  Figs.  4.3a-4.3d.  Also  from 
Figs.  4.4a-4.4d,  it  is  observed,  by  comparing  the  top  row  to  the  bottom  row,  that  the  strain 
energy  contour  lines  of  equal  value  all  move  towards  the  center,  and  that  the  contours 
shapes  of  the  strain  energy  over  GaAs  (001)  are  sharply  different  to  those  over  GaAs 
(111). 


Table  4.1  Maximum  strain  energy  £’max  on  the  surface  of  the  substrate  GaAs  for  different 
QD  shapes  with  different  depths  (unit  of  energy=l  18.8x1 


015Nm). 


Depth 

Cubic 

Truncated  pyramid 

Pyramid 

d=lnm 

(001) 

F 

^max 

2.54 

2.93 

2.71 

(*,y) 

(±1,±1) 

(±1,±1) 

(0,0) 

(111) 

F 

^max 

2.55 

2.75 

2.05 

0 x,y) 

(1,0) 

(1,0) 

(0,0) 

d=2nm 

(001) 

Emax 

1.77 

1.99 

0.96 

Cry) 

(0,0) 

(0,0) 

(0,0) 

(111) 

F 

^max 

1.31 

1.34 

0.75 

Cry) 

(0,0) 

(0,0) 

(0,0) 

d=3nm 

(001) 

F 

^max 

1.03 

1.09 

0.41 

Cry) 

(0,0) 

(0,0) 

(0,0) 

(111) 

F 

■^max 

0.70 

0.70 

0.36 

Cry) 

(0,0) 

(0,0) 

(0,0) 
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Figure  4.4  Normalized  strain  energy  on  the  surface  of  the  half-space  substrate  of 
GaAs  (001)  (top  row)  and  GaAs  (1 1 1)  (bottom  row)  induced  by  (a)  a  cubic,  (b)  a 
truncated  pyramid,  (c)  a  pyramid,  and  (d)  a  point  QD,  where  the  heavy  blue,  red 
and  black  lines  correspond  to  normalized  strain  energy  of  0.3,  0.6  and  0.9 
respectively.  The  QD  is  embedded  within  the  substrate  with  its  top  side  at  a  depth 
d=2nm  from  the  surface. 

We  now  study  the  effect  of  QD  depth  on  the  strain  energy  distribution  on  the  surface. 
The  top  row  of  Fig.  4.5  shows  the  depths  of  cubic  QD  within  the  substrate  where  d=lnm, 
2nm,  and  3nm,  whilst  the  strain  energy  distributions  on  the  surface  of  GaAs  (001)  and 
(111)  are  shown,  respectively,  in  the  middle  and  bottom  rows.  Again,  the  heavy  blue,  red 
and  black  contour  lines  correspond,  respectively,  to  the  strain  energy  values  of  0.3,  0.6, 
and  0.9.  It  is  apparent  from  Fig.  4.5  that  as  the  QD  moves  away  from  the  surface,  the 
strain  energy  contours  approach  to  those  due  to  a  point  QD  with  equal  volume. 
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Figure  4.5  Geometry  of  a  cubic  QD  within  a  half-space  substrate  with  different 
depth  d :  the  vertical  x-z  plan  view  (top  row),  and  the  strain  energy  induced  by  the 
cubic  QD  in  substrate  GaAs(OOl)  (middle  row)  and  GaAs  (111)  (bottom  row):  (a) 
d=  1,  (b)  d=2,  and  (c)  d= 3.  The  heavy  blue,  red  and  black  lines  correspond  to 
normalized  strain  enerev  of  0.3.  0.6  and  0.9  resnectivelv. 


Similar  numerical  results  are  shown  in  Figs.  4.6  and  4.7,  respectively,  for  the  truncated 
pyramidal  and  pyramidal  QDs.  While  the  strain  energy  distributions  induced  by  the 
pyramidal  QD  are  similar  to  those  by  an  equivalent-volume  QD,  those  by  cubic  and 
truncated  pyramidal  QD  are  different.  These  again  demonstrate  the  effect  of  the  QD 
shape  and  depth  on  the  strain  energy. 
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Figure  4.6  Geometry  of  a  truncated  pyramid  QD  within  a  half-space  substrate  with 
different  depth  d :  the  vertical  x-z  plan  view  (top  row),  and  the  strain  energy 
induced  by  the  truncated  pyramid  QD  in  substrate  GaAs(OOl)  (middle  row)  and 
GaAs  (111)  (bottom  row):  (a)  d=  1,  (b)  d= 2,  and  (c)  d=  3.  The  heavy  blue,  red  and 
black  lines  correspond  to  normalized  strain  energy  of  0.3,  0.6  and  0.9  respectively. 
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Figure  4.7  Geometry  of  a  pyramid  QD  within  a  half-space  substrate  with  different 
depth  d:  the  vertical  x-z  plan  view  (top  row),  and  the  strain  energy  induced  by  the 
pyramid  QD  in  substrate  GaAs(OOl)  (middle  row)  and  GaAs  (111)  (bottom row):  (a) 
d=  1,  (b)  d=2,  and  (c)  d= 3.  The  heavy  blue,  red  and  black  lines  correspond  to 
normalized  strain  energy  of  0.3,  0.6  and  0.9  respectively. 


4.5  Conclusions 

In  this  chapter,  we  presented  an  analytical  method  for  calculating  the  QD-induced 
strain  fields  in  half-space  semiconductor  substrates.  The  QD  is  assumed  to  be  of  any 
polyhedral  shape  which  can  be  efficiently  approximated  by  a  number  of  flat  triangles.  We 
studied  cubic,  pyramidal,  truncated  pyramidal  and  point  QDs  within  the  GaAs  (001)  and 
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GaAs  (111)  half-space  substrates.  The  numerical  results  illustrate  that  the  shape  of  the 
QDs  has  apparent  influence  on  the  strain  energy  distribution  on  the  surface,  so  is  the 
depth  of  the  QDs.  These  results  should  be  interesting  to  the  overgrowth  of  QDs  on  the 
substrate  where  the  long-range  strain  energy  on  the  surface  plays  an  important  role  in 
controlling  the  new  QD  patterns  and  shapes. 
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CHAPTER  V 


THREE-DIMENSIONAL  STRAIN  AND  PIEZOELECTRIC  FIELDS  IN  NITRIDE- 

BASED  QUANTUM  DOTS 


5.1  Summary 

Group  III  nitride -based  quantum  dots  structures  have  been  of  major  interest  in  recent 
years  because  of  their  ability  to  produce  high-quality  light  emitting  and  laser  diodes 
[ConOO,  Med02,  Not95].  This  has  been  largely  stimulated  by  the  progress  in  self- 
assembly  quantum  dot  growth.  In  order  to  achieve  uniform  dots  in  size,  shape  and  density, 
much  work  still  needs  to  be  done  to  understand  the  growth  of  these  structures. 
Theoretical  study  of  the  elastic  and  electric  fields  in  and  around  the  quantum  dots  is  now 
of  the  focus  interest  because  it  can  serve  as  a  useful  tool  to  guide  the  QDs’  growth  and 
also  to  understand  their  properties  [Jon06,  MakOl,  Zha05,  Pan07,  Pan08]. 

In  this  chapter,  a  simple  and  elegant  method  for  the  detailed  analysis  of  the  strain  and 
piezoelectric  field  in  electronic  structure  of  AIN  QDs  is  employed.  It’s  based  on  the 
Green’s  function  in  terms  of  Stroh  formalism  with  the  corresponding  exact  integration  of 
the  Green’s  function  over  the  QD  surface.  This  method  takes  into  account  quantum  dots 
of  arbitrary  shape  with  the  most  general  class  of  anisotropy  and  piezoelectricity. 

This  chapter  is  organized  as  follows.  In  section  5.2,  we  develop  the  method  to 
calculate  the  3D  fields  distributions  in  and  around  single  or  multiple  QDs.  The  results  are 
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given  in  section  5.3.  We  first  present  the  field  distributions  in  single  QD.  The  results  are 
in  good  agreement  with  previous  work.  After  that,  a  pair  of  QDs  are  studied.  Conclusions 
are  drawn  in  section  5.4. 


5.2  Problem  Description 

Assume  that  the  eigenstrain  y*  is  uniform  within  the  QD  domain  V,  the  induced 

displacement  at  point  d  can  be  expressed  as  boundary  integration 

uK(d)  =  CiJLmfLm  J  U j  (x; d )nt (x)dS (x)  (5.1) 

av 

Noted  that  Eq.(5.1)  is  similar  as  Eq.(4.1).  The  difference  is  that  we  use  the  extended 
notation  (which  is  described  in  chapter  II)  to  study  both  elastic  and  electric  fields. 
Uf(x',d)  is  the  /-th  Green’s  extended  displacement  and  can  be  expressed  as 

U(x;d)  =  U°°(x;d)+-^-Yf  AG,ATd6  (5.2) 

2  7T  ^ 

with 


-p,x3  +  Pjd3  -  [(.r,  -  dx )  cos  6  +  (x2  -  d2 )  sin  0 ] 


Follow  the  same  instruction  in  chapter  IV,  we  approximate  the  surface  of  QD  with  a 
number  of  flat  triangles  and  carry  out  the  integration  of  the  half-space  Green’s  function 
over  these  triangles. 

The  integration  of  the  complementary  part  over  a  flat  triangle  A  can  be  expressed  as 

UK  id)  =  C^rln.ji^  I’  AG,Ard#m(x)  (5.4) 
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In  order  to  find  the  extended  strain,  we  take  the  derivatives  of  Eq.  (5.1)  with  respect  to 


the  observation  point  d  to  obtain  the  QD  induced  strain  field 

Ykp (d)  =  \  rlmCuun  J  \U jJp  (x; d)  +  UPJA  (x;d)\n, ( x)dS (x) 
1  ay 


and  electric  field 

EP(d )  =  fLmCuLm  J  Ujdp(x;d)ni(x)dS(x) 

dV 


(5.5) 


(5.6) 


5.3  Results  and  Discussion 

We  now  apply  the  above  analytical  solutions  to  investigate  the  elastic  and  electric 
fields  induced  by  a  buried  QD  within  the  AIN  half-space  substrate.  We  consider  two 
different  growth  orientations:  one  is  AIN  (0001)  growing  along  the  (0001)  axis  and  the 
other  is  AIN  (1000)  growing  along  the  polar  direction.  The  material  parameters  used  are 
taken  from  [JogOl].  The  boundary  condition  is  assumed  to  be  traction  free  on  surface  and 
the  misfit  strains  are  yxx  =  y[y  =0.1375,  y  =0.1267  in  AIN  (0001),  and  yxx  =0.1267  , 

yyy  =  y  =0.1375  in  AIN  (1000).  A  schematic  3D  view  of  an  AIN  QD  is  shown  in  Fig. 

5.1(a)  with  its  cross  section  in  x-y  plane  shown  in  Fig.  5.1(b).  The  QD  size  is  determined 
by  the  top  and  base  diameters  of  the  pyramid,  i?b=8.5nm  and  i?t=4nm.  h=4. 1  run  is  the  QD 
height  and  d= 6.9nm  is  the  distance  between  top  of  QD  and  half  space  surface.  The  size  of 
the  QD  is  taken  from  experimental  values  [Arl99]. 
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(a) 

Surface 
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x 


Figure  5.1  Schematic  diagrams  of  hexagonal  truncated-pyramidal  QD  showing  dot  shape 
and  geometrical  parameters,  (a)  3D  view  of  a  single  QD  buried  under  half  space  substrate. 
Lines  A  and  B  parallel  to  z-axis  and  x-axis,  respectively,  denote  two  of  the  lines  along 
which  the  strain  field  has  been  plotted;  (b)  View  of  the  QD  structure  in  x-y  plane. 


Before  looking  at  multiple  dots,  a  single  isolated  hexagonal  truncated-pyramidal  dot  is 
studied  to  facilitate  comparison  with  previously  published  calculations  for  point  dot.  The 
strain  field  on  the  surface  is  an  important  part  of  the  overall  structure,  as  it  can  provide 
valuable  insights  into  the  understanding  of  the  growth  of  subsequent  layer  of  QDs. 
Shown  in  the  Fig.  5.2  are  the  contour  plots  of  hydrostatic  strain  on  the  surface  of 
substrates  AIN  (0001)  and  AIN  (1000)  due  to  a  buried  hexagonal  truncated-pyramidal 
QD  as  described  above.  From  Fig.  5.2,  one  can  notice  that  the  hydrostatic  strain 
distributions  are  symmetric  about  the  z-axis  and  the  surface  is  highly  strained.  Directly 
above  the  QD,  the  AIN  is  under  a  large  tensile  strain. 
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Figure  5.2  (a)  Contour  of  the  hydrostatic  strain  on  the  surface  of  AIN  (0001)  due  to  a 
hexagonal  truncated-pyramidal  QD.  (b)  Contour  of  the  hydrostatic  strain  on  the  surface 
of  AIN  (1000)  due  to  a  hexagonal  truncated-pyramidal  QD. 


Figs.  5.3(a)  and  (b)  show,  respectively,  the  contours  of  piezoelectric  potential  on  the 
surface  of  substrate  AIN  (0001)  and  (1000)  due  to  the  same  hexagonal  truncated- 
pyramidal  QD.  It  is  observed  that  the  contour  of  piezoelectric  potential  on  the  surface  of 
AIN  (0001)  is  rotationally  symmetric  about  the  z  axis,  and  exhibits  a  large  potential  well 
in  the  center,  which  will  contribute  to  the  lateral  carrier  confinement  in  the  dot,  while  the 
potential  on  the  surface  of  AIN  (1000)  is  symmetric  about  x  axis  but  anti-symmetric 
about  y  axis.  It  can  be  also  observed  that  the  magnitude  of  potential  on  the  surface  of  AIN 
(0001)  is  much  larger  than  that  on  the  surface  of  AIN  (1000),  i.e.  0.0065  VS.  0.0026. 
Those  results  in  Figs.  5.2  and  5.3  are  in  good  agreement  with  the  previous  published 
results  for  point  QD  [Pan03]. 
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Figure  5.3  (a)  Contours  of  the  piezoelectric  potential  on  the  surface  of  AIN  (0001)  due  to 
a  hexagonal  truncated-pyramidal  QD.  (b)  Contours  of  the  piezoelectric  potential  on  the 
surface  of  AIN  (1000)  due  to  a  hexagonal  truncated-pyramidal  QD. 


The  strain  field  yxx,  yyy  and  yzz  variation  along  line  A  of  Fig.  5.1(a)  are  shown  in  Fig. 
5.4(a).  From  Fig.  5.4(a),  it  is  seen  that  yxx=  yyy,  which  is  because  of  rotational  symmetry 
along  line  A.  Generally,  inside  the  QD,  the  forces  are  exerted  primarily  in  the  x-y  plane. 
As  the  substrate  tries  to  force  the  QD  lattice  constant  to  approach  that  of  substrate,  it 
induces  compressive  strains  in  the  x-y  plane  and  tensile  strains  along  z  axis.  The  signs  of 
the  strains  are  reversed  in  the  surrounding  substrate  because  of  the  opposing  forces. 
These  results  are  reflected  in  both  Figs.  5.4(a)  and  (b),  where  5.4(b)  shows  the  strain 
variation  along  the  surface  line  B  of  Fig.  5.1(b).  From  Fig.  5.4(b),  we  can  see  that  yxx,  yyy 
reach  maximum  directly  above  the  top  of  QD  which  makes  this  location  to  be  the  most 
favorite  for  the  incoming  adatoms  to  nucleate.  The  maximum  values  of  yxx  and  yyy  are 
equal  to  each  other  due  to  the  rotational  symmetry.  It  is  also  observed  that  yyy  does  not 
remain  tensile  along  the  entire  line.  It  changes  sign  and  becomes  compressive  away  from 
the  top  of  QD  due  to  energy  conservation. 
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z(nm)  x(nm) 

Figure  5.4  Line  scans  strain  distributions  in  AIN  (0001)  along  vertical  line  A  (a)  and 
surface  horizontal  line  B  (b)  in  Fig.  5.1. 

The  contours  of  piezoelectric  potential  in  the  x-z  plane  (y=0)  due  to  the  same 
hexagonal  truncated-pyramidal  QD  are  visualized  in  Fig.  5.5.  The  blue  (or  dark)  areas 
show  regions  of  low  potential  and  the  yellow  (or  light)  areas  regions  of  higher  potential. 
It  can  be  seen  that  there  is  a  large  potential  well  created  by  the  potential  difference 
between  the  pyramid  base  and  top  in  AIN  (0001)  while  large  potential  well  shows  along 
the  left  and  right  side  in  AIN  (1000).  The  vertical  profile  of  potential  in  AIN  (1000)  is 
anti-symmetric  about  z  axis.  We  can  also  see  that  the  potential  in  both  structures  are  in 
same  magnitude,  unlike  the  potential  on  the  surface  in  Fig.  5.3. 
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Figure  5.5  Contours  of  piezoelectric  potential  in  the  x-z  plane  (y=0)  in  (a)  AIN  (0001)  (b) 
AIN  (1000)  half-space  structures.  The  down  trapezoid  show  the  vertical  profile  of 

hexagonal  truncated-pyramidal  QD. 


Surface 


B 


Figure  5.6  Schematic  diagrams  of  two  coupled  hexagonal  truncated-pyramidal  QDs. 
Dx=25nm  is  the  distance  between  the  center  of  two  QDs. 


Fig.  5.7(a)  shows  strain  distribution  along  the  z-axis  shown  in  Fig.  5.6.  As  we  can  see, 
7xx  is  compressive  while  yyy  is  tensile  along  the  entire  line.  Fig.  5.7(b)  shows  strain 
distribution  along  the  x-axis  on  the  surface.  It  can  be  observed  again  that  yxx  does  not 
remain  tensile  along  the  whole  line.  It  is  also  noted  that,  unlike  the  single  dot  case,  the 
maximum  value  of  yxx  is  no  longer  equal  to  that  of  yyy,  which  means  the  rotational 
symmetry  for  a  single  dot  along  x-axis  does  not  hold  for  multiple  dots. 
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Figure  5.7  Line  scans  strain  distributions  in  AIN  (0001)  structures  due  to  coupled  QDs 
along  vertical  line  A  (a)  and  surface  horizontal  line  B  (b)  in  Fig. 5. 6. 


Figs.  5.8(a)  and  (b)  show,  respectively,  the  contours  of  piezoelectric  potential  on  the 


surface  of  substrate  AIN  (0001)  and  (1000)  due  to  the  pair  of  QDs.  It  is  observed  that  the 


contour  of  piezoelectric  potential  on  the  surface  of  AIN  (0001)  is  no  longer  rotationally 


symmetric  in  the  regions  directly  above  the  QD.  It  is  further  observed  from  Fig.  5.8(b) 


that  a  large  potential  well  is  created  between  the  two  QDs. 


x(nm)  x(nm) 


0.004 

0.0032 

0.0024 

0.0016 

0.0008 

0 

-0.0008 

-0.0016 

-0.0024 

-0.0032 

-0.004 


Figure  5.8  (a)  Contours  of  the  electric  potential  on  the  surface  of  AIN  (0001)  due  to  two 


hexagonal  truncated-pyramidal  QDs.  (b)  Contours  of  the  electric  potential  on  the  surface 


of  AIN  (1000)  due  to  two  hexagonal  truncated-pyramidal  QD. 
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5.4  Conclusion 


In  this  chapter,  a  detailed  theoretical  calculation  of  the  strain  and  piezoelectric  fields  in 
and  around  nitride -based  quantum  dots  buried  in  anisotropic  half-space  substrate  is 
presented.  The  method  is  applied  to  study  a  specific  example,  AIN  QDs  in  shape  of 
hexagonal  truncated-pyramidal.  First,  we  present  strain  and  electric  potential  distributions 
on  the  surface  induced  by  single  QD  and  compare  the  results  with  previous  published 
ones  for  point  QD.  Good  agreements  are  achieved.  The  vertical  profile  of  piezoelectric 
potential  in  single  QD,  as  well  as  the  strain  and  potential  in  a  pair  of  QDs  is  also  studied. 
We  find  that  the  pair  of  QDs  can  affect  the  symmetry  existed  in  isolated  single  QD  and 
also  the  large  piezoelectric  potential  leads  to  a  strong  lateral  confinement  in  the  QDs. 
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CHAPTER  VI 


CONCLUSIONS  AND  FUTURE  WORK 

6. 1  Conclusion  and  Contributions 

In  this  study,  the  elastic  and  electric  fields  in  or  around  QWR  embedded  in  anisotropic 
full-plane  and  QWR  free-standing/embedded  in  bimaterial-plane  are  investigated  by 
novel  Green’s  functions  and  BEM.  The  strain  energy  on  the  free  surface  of  half-space 
substrate  induced  by  cubic,  pyramidal,  truncated  pyramidal  or  point  QD  is  also  studied 
analytically.  The  main  conclusions  and  contributions  are  summarized  as  following: 

1.  An  accurate  BEM  modeling  is  presented  for  the  strain  and  electric  fields  analysis  in 

QWR  structures.  From  this  study,  some  important  features  are  observed. 

1)  Both  the  inclusion  and  inhomogeneity  models  predict  very  close  strain  field  in  the 
substrate  far  away  from  the  QWR. 

2)  For  points  inside  or  close  to  the  QWR,  the  strain  difference  between  the  inclusion 
and  inhomogeneity  models  can  be  as  high  as  10%,  which  can  result  in  strong 
variations  of  the  confined  electronic  states. 

3)  The  magnitude  of  strain  inside  the  QWR  is  much  larger  than  that  outside,  whilst 
the  electric  fields  have  the  same  magnitudes  both  inside  and  outside  the  QWR. 


83 


4)  While  the  difference  in  the  electric  fields  between  the  inclusion  and 
inhomogeneity  models  is  small  in  the  nitride  based  system  (InN/AIN),  that  in  the 
InAs/GaAs  (111)  can  be  large. 

5)  It  is  observed  that  the  electric  fields  in  the  QWR  depend  strongly  on  the  QWR 
geometry  shape.  The  electric  field  in  triangle  and  square  QWRs  is  sharply 
different  from  those  in  the  polygons  made  of  more  than  4  sides. 

2.  The  bimaterial  Green’s  functions  in  anisotropic  elastic  media  in  terms  of  the  Stroh 
formalism  are  derived.  It  is  then  applied  to  calculate  the  strain  energy  density  and  the 
relative  strain  energy  in  InAs  (111)  free-standing  QWR  or  in  QWR  embedded  in 
GaAs  (111)  substrate.  Several  interesting  features  are  found: 

1)  For  the  case  of  an  isosceles  triangle  of  InAs  QWR  on  the  substrate,  the  strain 
energy  density  within  the  QWR  strongly  depends  on  the  base  angle  of  the  triangle. 

2)  The  magnitude  of  the  relative  strain  energy  for  the  case  of  an  isosceles  triangle 
QWR  on  the  substrate  increases  with  increasing  base  angle. 

3)  For  a  square  QWR  either  free  standing  on  or  embedded  in  the  substrate,  the 
magnitude  of  the  relative  strain  energy  increases  with  increasing  depth  of  the 
QWR  in  the  substrate,  which  represents  the  competition  between  the  surface  and 
bulk  energies. 

3.  An  analytical  method  for  calculating  the  QD-induced  strain  fields  in  half-space 
semiconductor  substrates  is  developed.  The  QD  is  assumed  to  be  of  any  polyhedral 
shape  which  can  be  efficiently  approximated  by  a  number  of  flat  triangles.  We 
studied  cubic,  pyramidal,  truncated  pyramidal  and  point  QDs  within  the  GaAs  (001) 
and  GaAs  (111)  half-space  substrates.  It  can  be  concluded  that  the  shape  of  the  QDs 
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has  apparent  influence  on  the  strain  energy  distribution  on  the  surface,  so  is  the 
depth  of  the  QDs. 

4.  A  detailed  theoretical  calculation  of  the  strain  and  piezoelectric  fields  in  and  around 
nitride-based  quantum  dots  buried  in  anisotropic  half-space  substrate  is  presented. 
AIN  QDs  in  shape  of  hexagonal  truncated-pyramidal  is  studied.  We  find  that  the 
coupled  QDs  can  affect  the  symmetry  existed  in  isolated  single  QD  and  also  the 
large  piezoelectric  potential  leads  to  a  strong  lateral  confinement  in  the  QDs. 

These  results  should  be  extremely  useful  in  control  of  the  QWR  shape  and  size  in 
epitaxial  growth  and  simulating  QWR  semiconductor  structures,  and  in  the  overgrowth  of 
QDs  on  the  substrate  where  the  long-range  strain  energy  on  the  surface  plays  an 
important  role  in  controlling  the  new  QD  patterns  and  shapes. 

6.2  Recommendations  for  Future  Research 

Self-assembled  QDs  can  be  grown  layer  by  layer  to  form  ordered  nanostructures. 
Variations  in  sizes  and  shapes  of  self-assembled  QDs  have  been  a  critical  issue  for  device 
applications.  Multi-layering  of  self-assembled  QDs  can  provide  a  possible  route  for 
improving  the  size  uniformity  and  lateral  order  of  the  dots.  Theoretically  as  well  as 
experimentally,  in  multilayer  system,  the  elastic  interaction  of  dots  can  lead  to  the 
formation  of  long-range  spatial  correlations,  which  may  lead  to  a  lateral  ordering  and  size 
homogenization  of  the  dots.  This  self-organization  arises  from  the  preferred  nucleation  of 
the  islands  at  the  local  minima  of  the  strain  energy  on  the  surface  that  is  induced  by 
previously  buried  dots. 
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In  multi-layered  system,  different  types  of  vertical  correlations  have  been  observed. 
Wang  et.al.  [Wan04]  studied  the  vertical  correlation  (Figure  6.1a)  and  anticorrelation 
(Figure  6.1b)  as  a  function  of  the  spacer  layer  thickness,  where  vertical  anticorrelation 
means  that  the  subsequent  QDs  stack  in  the  space  between  the  underlying  QDs  and  but 
not  vertically  aligned  with  the  islands  below.  Springholz  et.al.  [SprOl]  treated  buried 
islands  as  point-like  defects  and  observed  the  ABAB... stacking  (Figure  6.2a)  and 
ABC  ABC...  stacking  (Figure  6.2b)  for  different  materials.  They  also  concluded  that  the 
type  of  vertical  stacking  crucially  depends  on  the  elastic  anisotropy  of  the  matrix  material 
as  well  as  on  the  growth  orientation. 


I - 1  I - 1 

40  nm  40  nm 


Figure  6.1  Cross-sectional  STM  images  of  the  Ino.5Gao.5As  QDs.  (a)  A  column  of 
correlated  QDs  in  the  75  ML  region  and  (b)  an  array  of  anticorrelated  QDs  in  the  150  ML 

region. 
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Figure  6.2  Schematic  illustration  of  the  possible  3D  dot  arrangement  in  self-organized  dot 
superlattices  (a)  Centered  tetragonal  arrangement  for  the  SiGe/Si  (100)  growth 
orientation,  (b)  trigonal  dot  arrangement  with  ABCABC:  :  :  dot  stacking  for  the 
PbSe/PbEuTe  (111)  growth  orientation. 

We  have  already  carried  out  part  of  this  research.  Using  our  exact  closed-form  Green’s 
function  and  taking  integral  over  all  the  QDs  surface,  we  calculated  the  strain  energy 
density  on  the  half-space  surface  induce  by  2x2  QD  arrays  in  InAs/GaAs  (001)  or  (111) 
system.  The  QD  is  assumed  to  be  cubic  and  has  a  size  of  10x10x5  (Figure  6.3). 
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surface 


Figure  6.3  2  by  2  QD  arrays  in  (a)  x-z  plane  and  (b)  x-y  plane 

Figure  6.4  shows  the  strain  energy  density  on  the  surface  in  both  (001)  and  (111) 
systems.  The  inter-distance  between  QDs  is  d=10.  The  red  squares  show  the  location  the 
buried  QDs  when  looking  topdown.  We  can  observe  that  when  the  depth  D  is  equal  to  2, 
the  preferred  nucleation  of  the  islands  will  be  located  exactly  above  the  buried  islands,  as 
shown  in  Figures  6.4a  and  6.4b;  when  the  depth  D  is  increased  to  20,  the  preferred 
location  for  new  QDs  will  be  at  the  space  between  the  buried  islands,  which  means  there 
would  be  a  vertical  anticorrelation,  as  shown  in  Figures  6.4c-6.4d. 
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Figure  6.4  Strain  energy  density  contour  on  the  free  surface,  a)  (001)  direction,  D=2;  b) 
(111)  direction,  D=20;  c)  (001)  direction,  D=2;  d)  (111)  direction,  D=20.The  distance 

between  QDs  is  d=10. 


If  we  the  inter  distance  d  between  the  QDs  is  decreased  to  5,  then  when  the  depth  of 
the  QDs  D=15,  we  can  also  observe  the  vertical  anticorrelation  (Table  6.1).  It  can  thus  be 
concluded  that  the  vertical  correlation  and  anticorrelation  are  function  of  the  spacer  layer 
thickness. 
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It  is  hoped  that  more  calculation  will  be  carried  out  on  the  elastic  and  electric  fields 
induced  by  larger  arrays  QDs,  such  as  4x4  and  10x10  arrays  within  the  multilayered  QD 
systems. 

Table  6.1  The  QD  depth  D  and  inter  distance  between  QDs  when  forming  vertical 

anticorrelation 


d  (inter  distance  between  QDs) 

5 

10 

D  (QD  depth) 

15 

20 
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